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ABSTRACT 

We thoroughly discuss a new supersymmetric grand unified sce- 
nario of the triplet seesaw mechanism where the exchange of heavy 
SU(2)w triplet states generates both neutrino masses and soft su- 
persymmetry breaking terms. This framework, recently proposed 
by us in a previous work, is highly predictive since it contains only 
three free parameters connecting low-energy neutrino parameters, 
lepton and quark flavour violation, sparticle and Higgs boson spec- 
tra and electroweak symmetry breakdown. These three parameters 
are the triplet mass Mr, the effective supersymmetry breaking scale 
Bt and a coupling constant A. We perform a complete analysis of 
the parameter space taking into account the present experimental 
constraints and considering different types of neutrino spectrum. 
A special emphasis is given to the particular features of the spar- 
ticle and Higgs spectra and to the model independent predictions 
obtained for the processes \i — > eX, fi — ► e conversion in nuclei, 
r — ► eY and r — > fiY (X = 7, ee, Y = 7, ee, fifi). In the appen- 
dices, we present some technical aspects relevant for our analysis. 
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1 Introduction 



It has long been recognised that the realm of neutrino physics may offer some insights on 
the search for physics beyond the Standard Model (SM). The evidence of non vanishing 
neutrino masses and of leptonic mixing angles, as provided by neutrino oscillations, calls 
for an extension of the SM particle content. Namely, the d — 5 operator Y„(LH) 2 /Ml, 
describing neutrino masses at the effective level [I], can be generated by decoupling some 
heavy degrees of freedom at the scale Ml where lepton number (L) is broken. This is the 
essence of the well-celebrated seesaw mechanism which can be realized at the tree-level by 
exchanging either singlet fermions N [2], or SU(2)w triplets T with zero [3] or non-zero 
hypercharge [4], at Ml- However, we should notice that, from a theoretical perspective, the 
concept of non-zero neutrino masses by itself does not tell us much about the new physics 
beyond Ml- On the other hand, alternative signals of lepton flavour violation (LFV) (besides 
neutrino oscillations) would be a clear and dramatic manifestation of new physics, since they 
are strongly suppressed within the SM by the smallness of neutrino masses. For this reason, 
and taking into account the increasing sensitivity of the present and future experiments, 
it is crucial to explore theoretical frameworks where such LFV processes can be sizeable. 
A typical example is the minimal supersymmetric standard model (MSSM) extended with 
renormalizable interactions that give rise to the d = 5 effective neutrino mass operator. 
Hence, at least two attractive features of the supersymmetric version of the above seesaw 
scenarios are worth to be recalled: 

1) supersymmetry (SUSY) alleviates the hierarchy problem of the SM, which would be 
exacerbated by the presence of one more high scale, Ml (besides the Planck mass scale 

M PL ) m 

2) lepton flavour violating (LFV) processes (otherwise unobservable) can be enhanced 
through one-loop exchange of the lepton superpartners if their masses are not too far from 
the electroweak scale and do not conserve flavour [6]. 

Regarding the latter aspect, most of the literature has been focussing on the most con- 
servative scenario of universal sfermion masses at a scale larger than Ml, which is realized 
in either minimal supergravity or gauge mediated supersymmetry breaking (GMSB) models 
with very large SUSY breaking mediation scale. In such cases, flavour non-conservation in 
the sfermion masses arises from renormalization group (RG) effects due to flavour-violating 
Yukawa couplings [7J [8j [9] encoded, at the effective level, in Y u . In this respect, it has been 
pointed out that, in the seesaw realization with non-zero hypercharge triplets, the flavour 
structure of the slepton mass matrix m| can be univocally determined in terms of the low- 
energy neutrino parameters [9]. In contrast, within the singlet seesaw, the determination of 
m| from low-energy observables requires model dependent assumptions [TOj, [TTJ [12] . 

Recently we have proposed a new supersymmetric scenario of the triplet seesaw mecha- 
nism in which the soft SUSY breaking (SSB) parameters of the MSSM are generated at the 
decoupling of the heavy triplets. Moreover, the mass scale of such SSB terms is fixed only by 
the triplet SSB bilinear term [13] • This scenario turns out to be highly predictive in the 
sense that it relates neutrino masses, LFV in the sfermion sector, sparticle and Higgs boson 
spectra and electroweak symmetry breaking (EWSB). In this work we aim to further eluci- 
date and discuss a more general version of this new framework including effects (previously 
neglected in Ref. [13]) which also entail quark flavour violation (QFV). 
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The paper is organised as follows. In Section [2] we review the main features of the SUSY 
triplet seesaw mechanism. Its embedding into the 577(5) gauge group is described in Section[3] 
where the role played by the heavy triplets as messengers of SUSY breaking is also discussed. 
It will be shown that the exchange at the quantum level of the triplet states generate all the 
MSSM SSB mass parameters. The results of the related analytical evaluations are presented 
and analysed in Section HI We proceed in Section with a detailed description of the 
flavour structure exhibited by the SSB terms. Then, in Section 0, we bring forward the 
phenomenological analysis of our framework. More specifically, we detail the experimental 
constraints imposed on the parameter space (Section l6.ip and present our numerical results 
(Section l6.2p . Section [7| is devoted to the relevant phenomenological predictions such as the 
sparticle and Higgs boson mass spectra (Section 17. ip and the expected size of lepton and 
quark flavour violation (Section 17. 2p . Afterwards, we discuss the predictions for several LFV 
processes (/i — > , r — > /ry , r — > ej , // — > eee , r — > jiee , r — > e^n , r — > r — > eee 
and fi — > e conversion in nuclei) and the peculiar correlations which arise in our scenario, 
together with a complete numerical analysis (Section I7.3p . Our conclusions and summary 
are drawn in Section [HJ Several technical aspects are collected in appendices: Appendix A 
is devoted to the generalization of the method based on the wave-function renormalization 
to derive the SSB mass parameters; Appendix B presents our analytical calculations of the 
MSSM coefficients for the LFV ijiiZ operators; finally, Appendix C regards the computation 
of the box diagrams relevant for the r — > /zee and r — > e/i/i amplitudes. 



2 Recalling the triplet seesaw mechanism 

Before starting the discussion of the main subject of our paper, let us briefly review the key 
features of the supersymmetric triplet seesaw mechanism. The requirement of a holomorphic 
superpotential implies introducing the triplets as super-multiplets T = (T°,T + ,T ++ ),T = 
(T°,f +,f ++) in a vector-like SU(2) w xU(l) Y representation, T ~ (3, l),f ~ (3, -1) 
The relevant superpotential terms are: 

-j^i%L i TLj + ^j=\ x H{TH x + -^=\ 2 H 2 TH 2 + M T TT + iiH 2 H u (1) 

where i,j = e,fj,,r are family indices, Lj are the SU{2) W lepton doublets and Hi(H 2 ) is the 
Higgs doublet with hypercharge Y = —1/2(1/2). The matrix Yip is a 3 x 3 symmetric matrix 
and Ai 5 2 are dimensionless unflavoured couplings. My and /i denote the mass parameters of 
the triplets and the Higgs doublets, respectively. 

By decoupling the triplet states at the scale Mr, one obtains the d — 5 effective operator 
Y U (LH 2 ) 2 /2M L where Y U /M L is identified as follows: 

1 Aq 

_!_Y« = ^-Y l i. (2) 

M L v M T T K ' 

At the electroweak scale the Majorana neutrino mass matrix emerges and is given bjQ: 

2 2 \ 

_ Mr, v ~ Mr T ' ( j 



1 The appearance of neutrino masses can be also interpreted in terms of non- vanishing vacuum expectation 
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where v 2 = v sin (3 = (H2) (v = 174 GeV). It is worth to emphasize that the flavour structure 
of the matrix Yy is the same as that of Y„ and hence of the neutrino mass . 

Without loss of generality, we choose to work in the basis where the charged- lept on 
Yukawa matrix Y e is diagonal. Therefore, all the information about low-energy lepton 
flavour violation is contained in Y„ or m,,: 



m,. 



mf = diag(mi,m 2 ,m 3 ) 



(4) 



where mi,m2,rri3 are the neutrino mass eigenvalues and the unitary lepton mixing matrix 
U can be written as 



U 
V 



V-diag(l,e^ 1 ,e i02 ) 



C12C13 

- Si2C23 — Ci2S23 s 13 e 
V s 12-523 — Ci2C23Si3e* 



iS 



C12C23 - 
— C 12"S23 



Sl 2 Ci 3 

— Sl2S23Si3e 
Sl2C23'Sl3e 



a 1 ) 



id 




(5) 



The mixing matrix V is responsible for LFV and, in particular, for neutrino oscillations. 
We have adopted the notation s^- = sinfly, c^- = cos^- for the three mixing angles 9 12 , 
623 and #i3, and denoted the "Dirac" and "Majorana" CP- violating phases by 5 and <f>i t 2, 
respectively. 

The relations ([3]) and (j4]) clearly show that the high-energy structure of Yy can be 
determined by the low-energy neutrino parameters (taking also into account the RG effects 
on the d = 5 operator which, however, do not introduce unknown flavour structures). The 
implications of such simple flavour structure become dramatic when one considers LFV 
induced by RG effects in the mass matrix m| of the left-handed sleptons [9j. Assuming, for 
instance, flavour- blind SSB terms at the grand-unification scale Mq (m| = m 2 ,!), the form 
of the LFV entries is 



m 



L> % 3 



tx mo(YlY T Wo£ 



M, 



G 



(6) 



which, in terms of the neutrino parameters, read 

2 



2 , M T \ f Mq 



A 2 ^2 



~ rrin 



M 



T 



M T ® \X 2 vi 



VfmfW 



loe 

H Mt 



(7) 



We note that the L-conserving combination YjY^ oc V(m^)Vt depends only on the neu- 
trino oscillations parameters, since the "Majorana" phases have been absorbed. From here, 
one can derive the relative size of LFV among the different leptonic family sectors [9] : 



(m|) TAI ^ 


v 




TfJ, 


(m|) re 


V(mf 


)2 V t" 


re 




[V 


; m ^) 2 vt 




( m l)^e 


[v(m?; 


2 V t] 





value (wew) induced by the EWSB on the scalar neutral state T , i.e. (T ) = -f^- In our scenario the values 
of (T°) will be smaller than 1CP 3 GeV, which is much below the upper bound of ~ 2 GeV inferred from the 
global fits of the electroweak data [TS] . 
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The above ratios depend only on the neutrino parameters which can be measured in low- 
energy experiments. This observation renders the SUSY triplet seesaw mechanism much 
more predictive when compared with the singlet one. Our scenario, therefore, constitutes a 
concrete and simple realization of the so-called minimal lepton flavour violation hypothesis 
which has been recently revived in the literature [TUd?]. In fact, relations like those shown 
in Eq. (jSJ) cannot be obtained in the latter case due to an ambiguity in the extraction of the 
high-energy neutrino parameters from low-energy observables [TO] . 

3 Gauge and Yukawa mediated SUSY breaking sce- 
nario 

A brief comment on the SSB pattern is now in order to motivate and introduce the main idea 
of this paper. We recall that the assumption of universality of the soft scalar masses at a 
high scale below Mpl may not be a justified one. As a matter of fact, flavour universality at 
Mpl, arising from some gravity- mediated SUSY breaking models [18] . is likely to be spoiled 
by Yukawa interactions in the energy range above Mq [19], [20]. Therefore, the common lore 
of adopting the universality condition at Mq should be regarded as a conservative approach 
to the issue of flavour violation (FV). In this work, we discuss an alternative scenario which, 
in our opinion, suggests a more motivated and predictive picture for the SSB pattern at high 
scales. This is the consequence of the fact that the triplet-exchange at the quantum level 
gives rise to the SSB term£| of the MSSM, i.e., the triplets play the role of SUSY breaking 
messengers [T3] . 

At first sight, the presence of extra SU(2)w triplet states at intermediate energies could 
spoil the simple gauge coupling unification, which can be achieved within the MSSM. How- 
ever, this can be remedied by invoking a grand unified theory (GUT) where the triplets 
live in a complete GUT representation, in such a way that gauge coupling unification 
can be preservecH. To this purpose, we arrange a SU(5) set-up where the T (T) states 
fit into the 15 (15) representation, 15 = S + T + Z with S, T and Z transforming as 
S ~ (6, 1, -§), T ~ (1, 3, 1), Z ~ (3, 2, |) under SU(3) x SU(2) W x U(1) Y (the 15 decom- 
position is obviousjl. 

The SUSY breaking mechanism is parameterized by a gauge singlet chiral supermultiplet 
X = Sx + dipx + 2 F Xl whose scalar Sx and auxiliary F x components are assumed to acquire 
a vev through some unspecified dynamics in the secluded sector (ipx is one of the goldstino 
components). In order to prevent the tree-level generation of SSB terms in the observable 
sector, the couplings of X to the MSSM fields must be forbidden. At the same time, X 
must couple with the 15 and 15 fields in order to trigger SUSY breaking in the messenger 
sector. Both these requirements can be achieved by, e.g., imposing that the SU(5) model 

2 In [21] the authors discussed the finite radiative contributions (arising from the decoupling of the triplet 
states) relevant to the SSB trilinear couplings in connection with the generation of the electric dipole mo- 
ments. 

3 This is not the only possibility to maintain gauge coupling unification in the presence of extra states at 
intermediate scales |22j . 

4 We find it interesting that such supersymmetric SU(5) with a 15,15 pair may be realized in contexts 
based on string inspired constructions [23l [24] . 
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Figure 1: The d = 5 supergraph from the 15, 15 exchange in the SU(5) symmetric phase (left). Once 
the coloured states t C 5# have been splitted from the Higgs doublet partners H2 and decoupled at 
the scale Mq, there is only the supergraph (right) with T,T exchange generating the L violating d = 5 
effective superpotential operator (LH2) 2 . 



conserves the combination of baryon and lepton number B — L = Q + ^Y, where Y are the 
hypercharges and 

Qw = g, Q5 = - g, Q5„(5„) = -g fg) , Q15 = g, Qte = g, = -2. (9) 

The SU(5) matter multiplets are understood as 5 = (d c ,L), 10 = (u c ,e c ,Q) and the Higgs 
doublets fit with their coloured partners t and i, like 5# = (t, H 2 ), 5h = {t, Hi). Given this, 
the relevant superpotential terms [13], consistent with the SU(5) and B — L symmetries, are 



1 

(10) 



W SU (5) = ^(Yi 5 5 15 5 + \5 H 15 5 H ) + Y 5 5 5 H 10 + Y 10 10 10 5 H + M 5 5 H 5 H + £X15 15 . 



The form of Wsu(5) makes explicit the fact that, thanks to the coupling with X, the 15 
and 15 states act as messengers of both B — L and SUSY breaking to the MSSM observable 
sector. Namely, while (Sx) only breaks B — L, (Fx) breaks both SUSY and B — L. These 
effects are tracked by the superpotential mass term M 15 15 15, where M15 = £{Sx), and the 
bilinear SSB term -B 15 M 15 15 15, with B 15 M 15 = -£(F X ). Once SU(5) is broken to the SM 
group@ we find, below the GUT scale M G [5], 

W = W + W T + W S , Z 

W = Y e e c H 1 L + Y d d c H 1 Q + Y u u c QH 2 + fiH 2 H 1 
W T = -^=(Y T LTL + \H2fH2) + M T Tf 

W s ,z = -^=Y s d c Sd c + Y z d c ZL + M Z ZZ + M S SS. (11) 
v2 

Here, Wq denotes the MSSM superpotential where the standard notation for the supermul- 
tiplets is understood. The terms relevant for neutrino mass generation [cf. (0Q)] are contained 
in Wt while the couplings and masses of the coloured fragments S and Z are included in 



5 For the sake of brevity, we have omitted in the SU(b) invariant superpotential (TTUj) other terms, as 
those involving the adjoint 24 representation responsible for the SU(5) breaking. 
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Figure 2: The d = 6 supergraph from the 15 exchange in the 577(5) symmetric phase (left). In the 
5£7(5) broken phase (/i < Mq), the latter splits into the remaining three supergraphes (right) which 
generate the B — L invariant d = 6 Kahler operators (LL)(ptf), (d c d c )(d c U^) and {d c L){d c] at 
low-energy. 

Ws y z- From the comparison of Eqs. (fTTj) and (pQ), we observe that the B — L invariance 
forbids the coupling XiHxTH% (which is not relevant for neutrino masses) leaving only the 
term proportional to A = A2. Consequently, the number of independent real parameters in 
Wt is reduced to eleven: Mx, A and the nine from the symmetric matrix Yy. 

Notice that we have relaxed the strict 577(5) symmetry relations for the Yukawa inter- 
actions and the mass term M5 by assuming S77(5) breaking effects (due to insertions of the 
24-representation), which are necessary to correct the relation Y e = Yj [25J and to solve the 
doublet-triplet splitting problem [26]. Thus, beneath the scale Mq, the coloured partners t 
and i are considered to be decoupled in order to adequately suppress dangerous d = 5 baryon 
number violating operators. The only d = 5 operator, generated by the exchange of the 15- 
fragments, is the L violating neutrino mass operator (see Fig. [1]). Instead, Fig. [2] shows that 
additional d = 6 operators arise from the exchange of T, 5* and Z [9 J . However, these are 
B — L conserving and, being suppressed by M T , are irrelevant for low-energy phenomenology. 
Hence, the presence of 5 and Z does not introduce new sources of baryon-number violation 
which could speed up the decay of the proton. 

As mentioned above, the B — L conserving superpotential of Eq. ( fTUl) contains the Ms- 
term from which the /i-parameter emerges in Wq. This mass parameter is not predicted 
by our model and will be determined through the requirement of EWSB (notice that the 
coupling XHiH 2 is forbidden by B — L). 

Also the coupling £ in Eq. (jTOj) could include 24-insertions, therefore allowing different 
masses for the 15-components. For simplicity, we consider the minimal case which implies a 
common mass M5 — — Mr — ^15 at the GUT scale. The Yukawa-unification condition 

Y S = Y T = Y Z = Y 15 , (12) 

can either hold or not at Mq, depending on the type and size of the 577(5) breaking effects. 
In this respect, we can discuss the following two extreme scenarios: 

(A) Eq. (fT2~|) holds at the GUT scale. As a consequence, flavour violation is extended to 
all the couplings Y5, Yt and Yz, originating FV effects both in the lepton and quark 
sector. In this work we shall focus on this general case. 
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(B) The states S and Z have negligible Yukawa couplings [Eq. f|T2|) does not hold because 
of SU(5) breaking]. In this case, FV is confined to the lepton sector since the couplings 
in Yt are the only ones which exhibit a non-trivial flavour structure. We have already 
studied this scenario in Ref. 1131. 



Beneath the messenger scale Mr, the particle content of our model is that of the MSSM 
and so, the superpotential reduces to the sum of the MSSM terms in Wq and the effective 
neutrino mass operator 

W <Mt = W + ^-(LE 2 )(LE 2 ). (13) 

As mentioned above, the vev of Fx induces the only tree-level SSB terms which, below the 
GUT scale, read (Bt = -B15 at Mq) 

- C>% T = B T M T {TT + ZZ + SS) + h.c. . (14) 

Such terms remove the mass degeneracy between the scalar and fermionic messenger com- 
ponents. To avoid tachyonic scalar messengers we require that £(Fx) < Mf> (or B? < Mr). 
At the tree-level, the ordinary MSSM supermultiplets are degenerate as they do not couple 
to the superfield X. Nevertheless, in the presence of the bilinear terms given in Eq. flHJ), 
the mass splitting is radiatively generated at the scale Mr through the gauge and Yukawa 
interactions between the messenger states S, T and Z and the ordinary MSSM fields. Our 
scenario can be, therefore, regarded as a gauge and Yukawa mediated SUSY breaking real- 
ization of the triplet seesaw mechanism. This will become clear in the next section where 
we present the complete MSSM SSB lagrangian. 

4 The SSB mass parameters 

We now discuss the SSB terms of the MSSM lagrangian £^f t SM which are generated at the 
decoupling scale of the heavy states S, T and Z. Below Mr, the most general £^jft SM can be 
written as: 

- £^f t SM =Ltm|Z + e c m| c e ct + &m 2 Q Q + dFm%d^ + u c m%u- ] + m\E\E x + m\ 2 E\E 2 
+ (E 1 e c A e L + E^AdQ + u c A u QE 2 + -M a X a X a + B H ^E 2 E 1 + h.c.) , (15) 

where we have adopted the standard notation for the slepton, squark and Higgs soft masses, 
the trilinear couplings A e rf u , the gaugino masses M a and the Higgs bilinear term B H . 

At the mass scale Mr, one-loop finite contributions are generated^ for A ej d jU , M a and 
B H . In Fig. [3] we draw the one-loop diagrams for A e , A d and A u (first and second rows), 
fiB H and M a (third row). The evaluation of these diagrams yields: 

A e = |^Y e (Y^Y r + Y^Y z ), 



167T 



6 Since the RGE- induced splitting on the masses of the 15-fragments is small, we decouple all the compo- 
nents at the common threshold Mt- 
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Figure 3: One-loop diagrams for the trilinear terms A e , and A u (first and second rows) and the 
bilinear Higgs term [iBh and gaugino masses M a (third row). The fields <3?(<i>) and ty($>) denote 
the scalar and fermionic components of the T(T), Z(Z), S(S) superfields in such a way that <I> = 

(T, Z, S% (T, Z), (5, Z) and § = (T, Z, S), (f, Z), (S, Z) for a = 1, 2, 3, respectively (the assignments 
for <l and ^ are obvious). 





T \ TT 



B T M T 



B T M T 




a„ - . . - a„ 

B t Mt - -* - B T M T 



vt 

Si -® 
B t Mt ' " ' B T M T 

z 



Figure 4: Examples of two-loop diagrams which contribute to m|. The upper diagrams generate 
contributions proportional to gf 2 (left) and to 9\2^t^t (right)- Instead, the lower diagram leads to a 
term proportional g^Y^Yz- 
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A« = SlA| 2 Y ?! 

l07T^ 

2B T 



167T 2 

16tt 

3i?r ,,|2 



Y Z Y^ + 2Y 5 Y+ )Y d 



^ = T^iAl 2 - (16) 

In these expressions, g a (a = 1, 2, 3) indicates the GUT normalized gauge coupling constants 
{91 — 92 — 93 at the unification scale) and the factor A^ = 7 in the r.h.s. of the gaugino masses 
M a corresponds to twice the Dynkin index of the 15-representation. A noteworthy feature of 
the above equations is that both the A-terms and the Higgs doublet bilinear parameter Bh 
are generated in our scenario^, once the Yukawa couplings and the \i parameter are present 
in the superpotentiaH-This is in contrast with the minimal GMSB realizations, where only 
the gaugino masses emerge at one- loop [23 EE]- We remark that the triplet states, being 
SU(3) singlets, would not communicate SUSY breaking to the gluino. This is one more 
motivation, besides the aforementioned gauge coupling unification requirement, to consider 
the GUT extension of the triplet seesaw with the 15 representation. 

Regarding the SSB squared scalar masses, the leading 0{F^/M^) = 0{B^) contributions 
vanish at the one-loop level due to a cancelation among the different terms [29]. For example, 
in the case of the soft masses m|, the one- loop diagrams driven by the exchange of the 
F-component of T cancel against those from the exchange of the F-component of L. Non- 
vanishing 0{Fj c /M^) contributions for the SSB squared scalar masses arise at the two- 
loop leve@. They are finite and can be evaluated either by diagrammatic computations 
or by means of a generalization of the wave function renormalization method proposed in 
Ref. [30j [31]. For the sake of illustration, we have depicted in Fig. 0] some representative 
two-loop diagrams which generate contributions to m|. The first diagram gives flavour-blind 
terms proportional to gf 2 - In turn, the remaining two diagrams generate LFV contributions 
proportional to g 2 2 Y^Yr and to gf 2 3Y z Yz- I n Appendix A we revisit the wave function 
renormalization approach and provide general formulas to extract all the SSB terms just by 
knowing the one-loop anomalous dimension matrices of the different fields above and below 
the messenger scale My. Our computation leads to the following result: 



mf 



\Bn 



21 * 21 4 /27 2 . nijWtv „ 2 _1_ O^ 2 _l_ 1 K„ 2 



Y^9l + y<?2 - I + 21a*J YJ,Y r - l-g( + 9gi + 16^ ) Y Z Y Z 



\1Qtt 2 / 

+ 18(Y+.Y T ) 2 + 15(Y^Y Z ) 2 + 3Tr(Y^Y r )Y^Y T + 12Yt,Y 5 Y+ Y z 
+3Tr(Y^Y z )Y^Y z + 9Y^Y|Y^Y t + 9(Y^Y r Y^Y z + h.c.) 



7 The SSB terms M a ,Bn and the trilinear couplings require interactions which violate the U(1)r sym- 
metry. In this context, the messenger bilinear terms (|14[) are those responsible for the U(1)r breaking. 

Although the \i parameter is not predicted by our model, we nevertheless assume that whatever mech- 
anism generates /i, it does not generate Bh- 

9 These 0{F\jM\) two-loop contributions dominate over the 0(F^/M%) = 0{B^jM\) one-loop ones 
for Mt > [AttYt / g 2 )Bj-- In the following such a hierarchy is fulfilled. 
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One can immediately recognize that the diagonal gauge-mediated contributions follow the 
expected form oc 2k a C[Ng^ (ki = 3/5, k 2 = k 3 = 1), where C[ is the quadratic Casimir of 
the /-particle [C{ = Y 2 = (Qf — T^) 2 and C = rL ^- for the fundamental representation of 
SU(n)]. 

The above expressions hold at the decoupling scale Mj- and are, therefore, meant as 
boundary conditions for the SSB parameters which then undergo (MSSM) RG running down 
to the low-energy scale fi s . All the soft masses have the same scaling property m so f t ~ 
Bt/{1Qtt 2 ) which, on the basis of the naturalness requirement m so f t ~ O(10 2 — 10 3 GeV), 
implies that Bt > 10 TeV. In principle, the SSB parameters also receive gravity mediated 
contributions of order the gravitino mass 1713/2 ~ F/Mpi (F 2 = {\Fx\ 2 ) + . . . is the sum of 
F-terms in the secluded sector). We assume such contributions to be negligible, which is the 
case if Mt -C 10 16 GeV £(Fx) / F. This also implies that the gravitino is lighter than the 
MSSM sparticles. 

It is apparent from Eqs. (flEI) and (jTTj) (and Figs. [3] and Hj) that the gauge interactions 
participate in the mediation of SUSY breaking in the gaugino and soft scalar masses. In the 
particular case of the sfermion mass matrices, the associated terms constitute the flavour 
blind contribution. Exactly the same situation occurs in pure GMSB models where flavour 
violation is automatically suppressed if the SUSY-breaking mediation scale is lower than 
the flavour or GUT scale [271 [2E1 [32]. In our case, SUSY breaking is also mediated by 
the Yukawa interactions Ys,Yt,Y.£ and A, giving rise to the bilinear parameter Bh, the 
trilinear couplings A e ^, u and to additional contributions to the Higgs scalar masses and the 
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sfermion mass matrices^- In this way, FV is transmitted to the SSB mass parameters. Since 
the origin of FV originally comes from the couplings Y15 (felt by the 5 supermultiplets), one 
expects that, due to the condition (j!2p . flavour violation is inherited with comparable size 
by the L and d c SSB parameters A e ,Ad, m| and rn| c . We emphasize that the FV entries of 
e.g., (m|)?j (i 7^ j) show up as finite radiative contributions induced by Bt at Mt, and they 
are not significantly modified by the running evolution to low-energy. This is different from 
a previous work [H] where a common SSB scalar mass m ~ (9(100 GeV) was assumed at 
Mq and the dominant LFV contributions to m| were generated by RG evolution from Mq 
down to the decoupling scale M T [see Eq. (JBJ)]. In such a case, finite contributions like those 
in Eqs. ([TBI) and (fT7|) also emerge at M T . Nevertheless, they are subleading with respect to 
the RG corrections, since Bt is of the same order as m . Instead, in the present picture, 
there is a hierarchy between the SSB parameter Bt and the remaining ones [see Eqs. ([TBI) 
and (HZ}], B\ > (B T g 2 /16n 2 ) 2 ~ m§ [T3]. 



5 Flavour structure in m| and m| c from neutrino pa- 
rameters 

In this section we discuss in detail the features of the flavour structure which emerges in the 
mass matrices m| and nr~ c , and in the trilinear couplings A e and A^. These parameters 
depend on the Yukawa couplings Y T , Y s , Y Z - The matrix Y T is determined at M T according 
to the matching expressed by Eq. (jSj) i.e., 

m„, 18 



2 



where the effective neutrino mass matrix m u has undergone the MSSM RG running from 
low-energy to Mt- The matrices Y5 and Yz are iteratively determined at Mt under the 
unification constraint of Eq. f[T2T) . For the determination of the mass matrix m u at low-energy 
[cf. Eq. OU)] we consider three different types of neutrino spectra: 

1. Normal hierarchy (NH): mi mi < 777,3, with 

m\ = m\ + Am|, 777.3 = m i + Am| + Am 2 A , (19) 



2. Inverted hierarchy (IH): 7773 mi < 7772, with 

777^ = TTlg — Am| + A777^, 7772 = m 3 + ^777^ , (20) 



3. Quasi degenerate (QD): 777,1 ~ 777 2 ~ 7773, with 

m 2 = m i + Am|, 7773 = 777^ + Ami + Am^, m\ ^> Am^ . (21) 
Other examples of Yukawa mediated SUSY breaking can be found e.g. in [29, ?33l [31j. 
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\m 2 — m|| and Am| 



The neutrino mass squared differences Am\ 
sponsible for the atmospheric and solar neutrino oscillations, respectively, together with the 
corresponding mixing angles 623 and 8u of the mixing matrix V [see Eq. fl5])]. From global 
analysis of the neutrino oscillation data, the following best fit values are available (with their 
allowed 3cr range) 



Am 2 A 
Am 2 s 



+1.1 

0.8 
1.0 
-0.8 



7.9 



10 
x 10 



-3 



eV 2 , sin 2 # 



eV 2 



23 : 

sin 2 On 



0.5 
= 0.3 



0.18 
0.16 
+0.1 
-0.6 



(22) 



For the lepton mixing angle #13 the following upper bound has been settled (at 3cr) 



sin 2 13 < 0.043 



(23) 



It is worth noticing that, for a NH spectrum, the largest entries of m„ (and, hence, of Yy) are 
those in the r — \i block, where the largest neutrino mass (Am^) 1 '' 2 dominate^]. Moreover, 
they are of comparable size due to maximal atmospheric mixing. On the other hand, the 
remaining matrix elements are smaller by one or two orders of magnitude, depending on the 
value of 6*i3 (which introduces the dependence on (Am^) 1 / 2 to these entries). In the IH case, 
while (m I/ ) ee and the r — // block are of the same order, (m^)^ and (m u ) Te are generically 
smaller. Finally, when the neutrino mass spectrum is of the QD type, the leading mass m x 
renders the matrix pattern less hierarchical and enhances the overall magnitude of m„ (or, 
equivalently, of Y T ). 

Let us now turn to the flavour structure of the SSB parameters which are generated at 
the messenger scale Mp. For simplicity of our discussion, we temporarily assume that the 
unification condition ffT2|) remains approximately valid at Mt, Yy ~ Y5 w Y#. In addition, 
we disregard the terms proportional to the Yukawa couplings Y e ^ in the expressions of m| 
and m 2 - c given in Eqs. ([171) . Under these assumptions, we can write: 



(A, 



(A rf )*- 



/ B 2 



V16vr 2 



3#r , v \ 

2 ^ [1 e)ii 



/ B 7 



V167T 2 



3B 



V(mf) 2 Vt 
V(m^ ) ) 2 V t " 
V(m?) 2 Vt 



r 



■in" 



ik 



ik 



(Yd) 



<9 2 J kj ~ 72 k 



c d a g 2 a 8 kj - 72 k 



V(m D ) 2 V t 



V(m^) 2 Vt 



6«Tr(m^) 2 <^ 



' M T 



6KTr(m?) 2 <^ 



(24) 



where ct , &„ are: c 



,i> c a are: c i _ ^ > c 2 _ 30 , C3 — 16, c\ — |§ , C2 — 6 and C3 — ™ 
of repeated indices is implicit). The flavour indices are i,j = e, /1, r (d, s,b) for the slepton 

> 2 and A„ 



m 



(squark) matrices m| and A e v , 
e <-> d, \i <-> s, r <-> 6 are understood. 



(summation 
• the slepton 

Wherever necessary, the correspondences 



n These considerations are valid for vanishing CP-phase, S = 0. From this point forward we restrict 
ourselves to this limit. The discussion about the impact of 6 ^ on our results is postponed to Section [7731 
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From Eqs. (jSJ)-©, arid considering the three different kinds of neutrino mass spectra 
shown in (fT9l)-fT2TT). we have: 



V(m?) 2 Vt 



mlSij + Am|V i2 V* 2 + (Am 2 A + Am 2 ) V l3 V* 3 (NH, QD) 
m 2 ^- + Am 2 A V t2 V* + (Am 2 A - Am 2 )V tl V* (IH) 



(25) 



which manifestly shows that the FV (i ^ j) entries are always independent of the lightest 
neutrino mass. The explicit expressions in terms of the neutrino parameters read: 



V(m?) 2 Vt; 
V(n0 2 Vt 
'V(mP) 2 V^ 



lie 



±Am\ C13C23S23 ± P (c 2 3S 2 3(Si2 - C 2 12 S 2 13 ) + C U S 12 S 13 {s 2 23 - C 

±Am\c l3 [s 13 s 23 ±pc 12 (c 23 si2 + C12S23S13)] 
±Am A c 13 [c 23 s 13 =F pc 12 (s 23 s 12 - C12C23S13)] , 



(26) 



where the upper (lower) sign applies for either the NH or QD (IH) neutrino spectrum, 
and p = Am 2 s /Am 2 A m 0.04. Due to the smallness of the parameter p, the quantities 
|(YyYr)y| oc |[V(m^) 2 V^]y| are not sensitive to the type of neutrino mass spectrum, 
unless S13 ~ PC12S12C23/C23 (or p C12S12C23/S23) ~ 0.02. For such values of s r3 , the pe and sd 
(re and bd) matrix elements in Eqs. (1241 would be strongly suppressed or even vanish for a 
IH (NH) mass pattern. Instead, the entries rp (bs) do not change significantly with # 13 , since 
the terms proportional to s 13 are suppressed by pcos(29 23 ) « 0. The quartic combinations 



V(m?) 2 Vt 



can be approximately obtained from the corresponding quadratic ones of 



)q. ( I26l) by performing the following replacements: 



NH : Am\ 
IH : Am\ 
QD : Am 2 A 



(Am 2 ) 2 (l + 2p) 



2 \2 



(Am Aj 
-> 2 m 2 Am 2 A 



p^p 2 (l-2p), 

P 
2 

P -> P- 



P 



2p 1 



(27) 



It is then straightforward to obtain the explicit expressions for the r.h.s. of Eqs. (1241) . taking 
into account Eqs. (126]) and ([271) . However, it is more instructive to consider the results for 
certain ranges of #13. Let us focus on m|, as the result for m 2 - c will follow directly. Distin- 
guishing the three types of neutrino mass neutrino spectra, we find: 

• NH spectrum, for s 13 <C p 2 C23Ci2Si 2 /s 2 3 ~ 8 x 10~ 4 : 



B 



T 



16tt 2 



B 



T 



16tt 2 



V(m^) 2 Vt 



V(mf?) 2 Vt 



tlx 



c e a g 2 a - 78 «A<(l + 2p) 
c^-6 K Ami(l + 14p) 



(2« 
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while for s 13 > PC23C12S12/S23 ~ 0.02: 



m 



2, ( v T v 



V(m^) 2 Vtl.. [c^ 2 -78«Am 2 (l + 2p) 



n 



IH spectrum for s 13 <C PC23C12S12/S23 ~ 0.02: 

(mi) T „ « -(^)^ [ V ( m ?) 2vt ] T , (^-84*^), 



while for s 13 > 2pc 23 Ci2Si2/s23 ~ 0.04: 



(29) 



(30) 



V167T 2 / 

QD spectrum (for any value of Si 3 ) 

(m ^"-(l^) 2/t [ V(m " )2Vt 



'J 



'^-162^(1 + ^ 



(31) 



(32) 



The above expressions exhibit the relevant flavour violating factors \V {m®) 2 V\j separated 
from the piece composed by terms proportional to the gauge couplings g\ and terms propor- 
tional to k Am\ (NH and IH) or nrnf (QD). For a given effective SUSY breaking scale B?, 
the overall size of those entries is controlled by k oc (M^/A) 2 . It is interesting to notice that, 
for a certain value of M T /X, the g 2 -terms are canceled and thus all the three off-diagonal 
entries can be vanishing. If S13 0.02, this peculiar flavour-blind suppression happens when 



M t/A ~ ^ 2 (t^) 1/2 and M T /\ 



~ v 



2 v 84Am 2 . 



I 1 / 2 , for the NH and IH neutrino spectrum, 



2( c%gl \i/2 



162mf 



, irre- 



respectively. Instead, for the QD case, the suppression occurs for M T /\ ~ f| 

spective of Si 3 . If Si 3 is tiny and the neutrino spectrum is either NH or IH, the entry (rp) is 
suppressed for Mr/A smaller (larger) by ps 1/3 (~ 1.4), with respect to the entry (pe) [(re)] 
for the NH (IH) spectrum. We can summarise this analysis saying that the off-diagonal en- 
tries of m| and m 2 ~ c can undergo a strong reduction in two cases: 1) when Si 3 ~ 0.02 and so 
the quantities [V(m^) 2 V^] Te and [V(m^') 2 V'l'] /je are vanishing for the NH and IH spectrum, 
respectively; 2) when there is a cancelation between the quadratic and quartic contributions, 
which depends on the parameters and A and on the neutrino spectrum. 

In the parameter range where the FV entries (m|)y are dominated by either the quadratic 
or the quartic Yukawa terms, the relative size of LFV (QFV) in the p — r (b — s) and p — e 
(s — d) sectors, does not depend of the ratio M T /\, and can be approximately predicted in 
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-^23/ 12 
-Rl3/12 



observablee 


3 V and 


(m|) T „ _ 


(m| e ) 6s 




(m| c ) sd 


(m|)re _ 




(m|) Me 





V(mf) 2 Vt 



Jr/i 



[V(m?)2Vt] Me ' 
V(m?) 2 Vt 



[V(m?) 2 Vl 



(33) 



Plugging the results of Eqs. fl26]) into the above expressions, and taking the central values 
( 1221) for the neutrino parameters, we have: 



#23/12^3=0 « 40 (-40), i? 2 
^13/12| si3=0 ~ — 1 ( — 1) #13/12 |si 3 =0.2 

for the NH (IH) spectrum. The QD case gives the same results as the NH one. 



3.2(3.8), 
0.8(1.2), 



(34) 



6 Phenomenological viability 

In the previous sections, we have set up a theoretical framework consisting of the super- 
potential (jl~3l) and the SSB lagrangian (fT5|) . beneath the energy scale M T , whose relevant 
soft mass parameters (|T6| [171) emerge at the messenger energy scale. Next, we describe the 
criteria employed to constrain the parameter space spanned by Mt, Bt and A. With the 
purpose of investigating the phenomenological viability of our scenario, a detailed numerical 
analysis is carried out in Section 16.21 The results are then thoroughly discussed. 



6.1 The phenomenological constraints 

Our approach to relate the low-energy measured parameters with the high-energy quantities, 
such as the Yukawa couplings, follows a bottom-up perspective. In particular, this allows us 
to determine the following quantities: 

• The matrix Y^ (as well as Yj and Y#) is determined at as already explained in 
the previous section [see Eq. ({TBI) ]. 

• The Yukawa matrices Y e , Y u , Y^ are extracted from the corresponding charged fermion 
masses, modulo tan (5. 

• The remaining low-energy input parameters /i(// s ) (with its sign) and tan/3 are deter- 
mined by the EWSB conditions: 

9 - tan 2 3m 2 Hn + m 2 N Ml n 2uB H , . 

2 = ^H2 Hi z sin 2/3 = —5 . 35 

P tan 2 /3-l 2 ' h m 2 Hl + m\ + 2fi 2 K ' 



12 We recall that, for low/moderate values of tan/3, the off-diagonal entries ([24)) are not substantially 
modified by the RG running from the high scale Mt to low-energies. Therefore, in Eq. ([33]) we have 
disregarded the renormalization effects on the neutrino mass [see Eq. ([15]) ]. However, this does not alter the 
form of m„, since it amounts to an overall correction factor, except possibly for the entries (m„) T i which 
receive extra (overall) corrections in the regime of large tan (3. 
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Table 1: Present limits and future sensitivities for the branching ratios (BR) of several LFV 
processes. The bound on the /i — > e conversion rate (CR) in Ti is also shown. 

We have included the one-loop tadpole corrections t\ t 2 through the redefinition of the 

Higgs scalar masses, i.e. rn 2 N = m\r 7=-^ — ^,m 2 Mn = m 2 Hn / - t2 . - . According 

to the standard practice, the above minimisation conditions are imposed at the scale 

Hence, we are left with three free parameters, Mt, Bt and the coupling A, which can 
all be taken to be real, without loss of generality. The parameter space spanned by these 
parameters is constrained as follows. 

• We impose the constraint coming from the decay fi — > ej, using the results in [TTJ. The 
experimental upper bound for the branching ratio of this decay is shown in Table [TJ 
The relevant LFV entry is (m|) Me [see Eqs. (IMl) and (126]) ]. 

• A conservative limit upon the lightest Higgs boson mass, nth > 110 GeV, is consid- 
ered on the basis of the negative LEP2 direct searches [17]. Our predictions on rrih 
include the low-energy radiative corrections which are obtained linking our code to 
FeynHiggs [48J. 

• We also require that the sparticle masses (which pass the tachyonic test) respect the 
experimental lower bounds set by Tevatron and LEP direct searches |15j . 



The SUSY contribution to the anomalous magnetic moment of the muon Sa^ [TT] has 
been considered and subjected to the following constraint at 99% C.L. 



- 13 x 10" iU < 5a^ < 57 x 10~ iU . (36) 

However, we will see that this requirement does not imply an important constraint on 
the parameter space. 
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• We have checked that all the coupling constants remain in the perturbative regime 
up to the scale Mq. In particular, while the presence of a complete GUT repre- 
sentation below Mq does not alter the value of Mq, the gauge coupling go gets an 
additional contribution bg^? = |^ln(M<r / 'Mq). The requirement of perturbativity im- 
plies ^ln(M G /M T ) < 24 which, for N = 7, gives M T > 10 7 GeV. The parameter 
space is further restricted by imposing the perturbativity requirement on the coupling 
constants A, Y T , Yg and Y Z - 

As it will be shown in the next section, the most stringent constraints come from the 
sparticle spectrum, fi — > cy decay, the Higgs boson mass and the requirement of radiative 
EWSB and perturbativity. 

6.2 Parameter space analysis 

In Fig. [5] the constraints imposed on the parameter space (A, My) are shown for = 
20 TeV and Sis = 0(0.2) in the upper (lower) panel, taking the NH neutrino spectrum 
(later we will comment on the other cases). The light-grey regions are excluded by the 
perturbativity requirement. For each value of Mt, there is a minimum value of A, scaling as 
~ 2 x 10~ 4 (Mr/10 11 GeV), below which the couplings Yy and/or Y reach the Landau 
pole between Mt and Mq. Similarly, there is an upper bound for A, above which A itself 
blows up. 

The EWSB constraint excludes the red region (which covers the range with A ~ 1 along 
the whole interval of Mt) limited on the left by the least achievable value of tan (3 ~ 2.5. We 
have studied the sensitivity of our results to the top mass, considering the la allowed range 
m t = (172.5 ± 2.3) GeV [50]. The most relevant effects of varying m t are the ones induced 
by the top- Yukawa term in the RG running of the soft mass J7Vg- . However, this variation 
does not affect our results considerably and, therefore, we will mainly take the central value 
m t = 172.5 GeV. 

The EWSB conditions select values of tan f3 (thick solid-contours) up to moderate one£l, 
(tan/3 < 25). Notice that the largest values of tan/3 are achievable only for M T < 10 9 GeV 
and small values of A < 0.2. In fact, Eq. (|35|) shows that larger values of tan/? (or smaller 
sin 2(3) require suppressed Bh{^ s ), which can be obtained by shortening the running energy 
interval (and so, smaller Mt are needed). We also display the iso-contours of the fi(^s) 
parameter (dashed-lines) as obtained by the EWSB conditions. One can realize that this 
parameter slightly increases with Mt due to the enhancement of the RG factor which affects 
m 2 H2 (fi s ) in the minimisation condition of Eq. (I3"51) . The whole allowed parameter space 
covers the range fi ~ (450 — 550) GeV. 

To better understand the behavior of the \i and tan (3 contours, it may be useful to con- 
sider Fig. [5] together with Fig. Indeed, the latter displays the solutions of Eq. (1551) for \i 
(upper plots) and Bh (lower plots), when Mt = 10 9 (10 13 ) GeV in the left (right) panel. We 
have chosen four distinct points in the (A, Mt) parameter space and shown, for each case, 
the predicted parameter Bh{^ s ) = Bh{Mt) + ABh, where Bh{Mt) is determined from 
Eqs. ( IT6|) and ABh englobes the RG running between Mt and fi s . The curves of Bh(i^ s ) 

13 As already mentioned above, varying m t in its 1 a interval does not lead to appreciable changes on the 
values of fj, and tan/3. Thus, we have only displayed the tan/3 and [i iso-contours for the central value of m t . 
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Figure 5: The (A, Mr)-parameter space configuration for Bt = 20 TeV and S13 = (S13 = 0.2) in 
the upper (lower) panel. Perturbativity and EWSB exclude the light-grey and red regions (see legend 
in the upper panel), respectively. The orange region is excluded by the Higgs mass bound < 110 
GeV for m t = 170.2 GeV, while the left-most (right-most) dash-dotted line delimit the same region 
for m t = 172.5(174.8) GeV. Inside the yellow and green areas, BR(/x — ► erf) is above the present 
experimental upper bound and the lightest slepton mass is below 100 GeV, respectively. In the dotted 
and hatched areas, the neutralino Xi iS the lightest MSSM sparticle and l\ is tachyonic, respectively. 
The thick-solid (dashed) lines correspond to the isocontours of tan/3 (// in GeV) for m t = 172.5 GeV. 
The thin solid lines refer to the FV parameter 5f s [defined in Eq. (j4§|) ]. See text for more details. 
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Figure 6: The same as in Fig. 5 for Bt = 50 TeV. 
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Figure 7: Left panel: Curves of Bh (dashed lines) obtained at fi s from the running of Bh(Mt), 
and of Bf, w (solid lines) determined by the minimization conditions ( 1351) , as a function of tan/3 
for B T = 20 TeV, M T = 10 9 GeV and A = 1(T 4 ,0.75 (lower plot). In the upper plot, the 
corresponding values of the fi parameter are shown. Right panel: The same as in the left one 

but for M T = 10 13 GeV and A = 0.25, 1.0. 



have to be compared with those of Bh as extracted from Eq. (|35|) (denoted by B^ in the 
plots). The crossing of these two curves (indicated by a black dot) signals the presence of 
a solution at the corresponding value of tan j3. For instance, the left-panel of Fig. [7] shows 
that for B T = 20 TeV, M T = 10 9 GeV and A = 10- 4 (0.75), EWSB occurs for tan/3 w 22(3) 
with /i w 450(400) GeV. A similar example is illustrated in the right-panel of Fig. [7] for 
B T = 20 TeV, M T = 10 13 GeV and A = 0.25 and 1. 

Consider now the constraint imposed by the lower bound on rrih- The orange region in 
Fig [5] shows the portion of the parameter space forbidden by the condition nih > 110 GeV, 
taking the lower limit of the la range for the top-mass (mt = 170.2 GeV). The depen- 

4 

dence on m t comes from the low-energy radiative corrections oc — ^\ti(u 2 I m%) . As m t in- 

creases, the region shrinks, being delimited by the left-most (right-most) dashed-line when 
m t = 172.5 (174.8) GeV. It is also clear that, for each value of M T , the upper bound on 
A is set by the Higgs mass constraint e.g., A < 0.45 — 0.55 when M T = 10 9 GeV. Since 
the dependence of the Higgs mass on tan/? mostly comes from the tree-level contribution 
| cos2/3\Mz, the iso-contours of m^ closely follow those of tan/3, setting in this way the least 
allowed value of tan (3 ~ 7. 

In our framework, the lightest MSSM supersymmetric particle is tipically the lightest 
slepton ii, except inside the dotted region where m^o < m^ (which is almost entirely ex- 
cluded by the Higgs and fi — > e'j constraints). The mass of l\ turns out to be below the LEP2 
lower bound of about 100 GeV in the region of the parameter space filled in green, where 
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the large values of tan (3 > 22 reduce rug through the left-right mixing at the electroweak- 
symmetry breaking (see also Section [7TT|) . In the hatched region (lower- left corner), where 
tan (3 ~ 30, the lightest slepton is tachyonic. 



All the constraints discussed so far, being related to 'unflavoured' observables, are not 
sensitive to the angle 6*13, as it is shown by the comparison between the upper (sin# 13 = 0) 
and lower (sin # 13 = 0.2) panels of Fig. On the contrary, the size of LFV strongly depends 
on it [see Eq. ( 126}) ] and hence, the region excluded by the bound on BR(/i — > ej) changes 
when different values of #13 are considered. For Mt ^2x 10 8 GeV, the /1 — > constraint 
provides the most restrictive lower bound on A. This stems from the fact that the size the 
LFV entry (m|) Me scales as (Mr/A) 2 [Eq. ff24l) ]. Consequently, the allowed A range is wider 
for lower values of Mt, closing up for Mt ~ 4 x 10 13 GeV. By switching on S13 (lower panel), 
the size of fie LFV is enhanced as 



and, correspondingly, the fi — > ey limit implies that the lower bound on A increases by a 
factor of ~ (140) 1//4 ~ 3.5 and so, the allowed parameter space shrinks. Taking into account 
all the constraints considered above, one concludes that values of Mt ^2 x 10 13 GeV are 
excluded for B T = 20 TeV, independently of the value of A. 

An equivalent analysis is presented in Fig. [6] for Bt = 50 TeV. The comparison with 
the previous case indicates that the regions excluded by the perturbativity and EWSB re- 
quirements are not significantly affected. On the other hand, the constraint on the Higgs 
mass implies a slightly different upper bound on A (A < 0.6) and a decrease of the mini- 
mum allowed value of tan/3 (tan/3 > 4). Indeed, the sparticle spectrum is now heavier and 

2 

thus, the radiative corrections ~ m (^0 to are larger. For this reason, smaller values 
of tan/3 are tolerated in the tree- level contribution ~ Mz\ cos 2(3 |. Contrarily to the previ- 
ous case, l\ is never tachyonic and its mass lies above the LEP bound. Consequently, for 
Mt = 10 7 — 10 s GeV the allowed A-range is much more extended with respect to the case 
with B T = 20 TeV. Moreover, larger values of tan/3 are now possible (tan/5 < 40). The 
heavier spectrum also makes the fi — > e^y constraint weaker, reducing the excluded yellow 
region. In conclusion, the allowed parameter space enlarges when Bt increases. Also for this 
case, the effect of non-zero #13 (lower panel) produces similar results as for smaller Bt, by 
raising the lower bound of A. 

Before concluding this section, a comment is in order about the influence of the type 
of neutrino spectrum on the allowed parameter space. In the IH case, the perturbativity 
constraint on the Yukawa couplings Yt,s,z would just imply a slightly larger minimum of 
A, for each Mt- All other constraints would instead be unaffected. For a QD spectrum, 
the effect from the perturbativity requirement would be much stronger (depending on the 
magnitude of the overall neutrino mass mi) since, as already mentioned, all the Yt,s,z entries 
increase when compared to the NH or IH cases. Therefore, the light-grey region would mostly 
cover the yellow area excluded by the BR(/i — > ej) bound. In conclusion, either for the IH 
or QD spectrum, the resulting allowed parameter space would not be much different from 
the NH case displayed in Figs. [5] and [6j 




(37) 
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7 Phenomenological predictions 

After describing the main phenomenological constraints imposed on the (A, My) parame- 
ters space, we will now go through the specific features of the sparticle and Higgs spectra 
(Section 17. ip of our scenario, and to the implications for several low-energy LFV processes 
(Section O). 



7.1 Sparticle and Higgs spectroscopy 

The spectrum of the superpartners is determined by the finite radiative contributions to 
the SSB parameters at M T (see Section HJ), acting as boundary conditions, and by the 
subsequent MSSM RG running from M T to /x s . The physical scalar masses are obtained by 
taking into account the latter effect at one-loop level and the low-energy D and F-term 
contributions. To get some intuition on the main features of the physical spectrum, we 
present some qualitative arguments in addition to the complete numerical results. 

In the present framework, the boundary conditions for the gaugino and sfermion masses 
are not universal at Mr, even when Mp is not far from Mq. This is due to the different 
gauge quantum numbers of the MSSM field representations [see Eq. (fTTl) ]. At lowest order, 
the gaugino masses M a at the messenger scald 1 " 4 ! are in proportion to the gauge coupling 
squared, M\ : M 2 : M3 = a\ : 0.2 '■ «3 {oi a = g\jAm). This relation is maintained at 
low-energy, like in unified SUGRA scenarios. The low-energy gaugino masses are given by: 

MM = ^Ng 2 a (ii s ), a = 1,2, 3, (38) 

which leads to: M 2 « 1.86 Mi « 0.35 M 3 (taking fi s ~ 700 GeV). The most interesting 
aspect comes from the fact that the scalar and gaugino masses are mutually related at 
the messenger scale. For the sake of our discussion, let us disregard for the moment the 
contributions proportional to the Yukawa couplings in the expressions of the scalar masses 
in Eqs. (iTTj) . as well as the Yukawa effects in the (1-loop) renormalization. In such a case, 
the low-energy sfermion masses have the form: 

where b\ = 33/5, b 2 = 1,63 = —3. The first term in nr~(/i s ) corresponds to the high-energy 
boundary contribution while the second one, AR a , accounts for the RG effects induced by 
the gaugino masses. The squark masses mq receive the main RG correction from the gluino 
mass term, which amounts to a positive shift on mj~. The term Ai? 3 is larger than the 
boundary condition value (since |1 — x%\ > \b a \/N = 3/7), by a factor of approximately 
8 (2) for M T = 10 14 (10 7 ) GeV. Notice that, in the minimal GMSB model (N = 1), the 



'in this section, we denote by overbar any quantity evaluated at the scale Mt- 
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Figure 8: The sparticle and Higgs spectrum for Bt = 20 TeV as a function of My and correspondingly 
taking A (as indicated in the upper horizonal axis) along the isocontour BR(/i — ► cy) = 1.2 x 10 -11 . 
Left panel: physical squark masses, m„ (black solid lines), mj (red dashed) and the pole gluino mass 
m~g (blue dash-dotted). Right panel: charged slepton masses (solid black lines), sneutrino masses 
my (dotted blue), neutralino and chargino masses m^o and (dashed blue). The dash-dotted green 
lines correspond to the Higgs boson masses m/j,m#,TOA and m H ±. In the small inner plot (on the left 
panel) we show the behavior of tan/3 (red line) and \x as obtained by the EWSB conditions. 



dominance of AR3 holds only for a messenger scale above 
expect the low-energy ratio ttlq/M^ to be given as 



M 3 



2Cf 1 - x 3 



1 + 



N{l-x 3 ) 



10 u GeV 



1/2 



Therefore, we 



(40) 



10 



11 



10 7 GeV. Hence, the gluino is the heaviest of 



which lies within 0.9 — 0.8 for 
the coloured sparticles. 

In Fig. [S] we display the sparticle and Higgs spectra for B T = 20 TeV. The values of A 
(shown in the upper horizontal axis) and Mr , are the ones along the border which delimits 
on the right the fi — > cy excluded regioro (see upper panel of Fig. [5]). The gluino pole 
mass includes the low-energy finite corrections (see e.g. [53]), which ammount to 20 — 30% 
of the tree-level value M3. As anticipated the gluino is the heaviest superparticle. The first 
and second generation squarks d (dashed curves) and u (solid) (mainly composed by left- 
handed squarks) are the heaviest. Instead, the lightest squark t\ is mainly a right-handed t c 
whose mass is pushed down by a negative shift, driven mainly by the top Yukawa-induced 
renormalization. This effect, together with the left-right squark mixing, is enhanced for 
Mt ^ 10 13 GeV where tan f3 becomes smaller (see the inner panel where we have drawn the 



15 The sparticle and Higgs boson spectra do not depend much on A and so, the results shown in Fig. [8] are 
quite representative of the whole parameter space at Bt = 20 TeV. 
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behavior of tan/3 and /i versus Mt). This analysis has shown that for Bt = 20 TeV in 
the allowed (A, Mt) portion, where the \i — > is close to the present bound, the squark 
spectrum lies in the range 800 — 950 GeV. This mass range will be soon explored by the 
Large Hadron Collider (LHC) [51] with a luminosity £ = 100 pb™ 1 and a center-of-mass 
energy y^s = 14 TeV. 

The left-handed slepton masses can instead be compared with the SU{2) W gaugino mass 
M 2 . In this case the renormalization factor in Eq. (1591 is comparable to the boundary 
contribution. At low-energy, the ratio m^jM^ lies in the range 1 — 0.5 for Mt = 10 14 — 
10 7 GeV. In the right panel of Fig. [8]we show the physical spectrum of the electroweak states. 
The heaviest charged sleptons £4,5,6 (which are mainly left-handed sleptons) have masses 
around 250 — 320 GeV in the allowed range Mt ~ 10 8 — 10 13 GeV. The sneutrino masses 
are splitted from the latter states mainly by the SU(2) W D— term, m| — m~ = cos2/3|. 

The right-handed slepton masses receive contributions only from the U(l)y gauge in- 
teractions and, therefore, are smaller than the other sfermion masses. In the range Mt = 
10 14 — 10 7 GeV, the high energy contribution is larger than the renormalization term AR% 
by a factor of ~ 1 — 3. The low-energy ratio m~ c jM\ is approximately in the range 1 — 0.7 
for Mt = 10 14 — 10 7 GeV. To these estimates one has to add the negative shift from the Y T - 
induced renormalization and the effect from the left-right slepton mixing at the electroweak- 
symmetry breaking. Both these contributions are important for large tan f3 and lower the 
physical fi mass below rrifc. In Fig. [8] (right panel) we can see that for Mt ^ 10 8 GeV the 
f\ mass is pushed below 100 GeV and, in the range Mt ~ 10 s — 10 13 GeV, f\ is indeed 
the lightest MSSM sparticle. For larger values of Mt, corresponding to small tan/?, the 
neutralino X\ becomes the lightest MSSM sparticle. Hence, either fx or x\ would decay into 
the gravitino, which is in fact the lightest SUSY particle in our framework. In conclusion, 
the slepton masses lie in the range ~ 100 — 320 GeV (for Bt = 20 TeV) and, therefore, are 
within the discovery potential of the LHC. 

Concerning the physical charginos and neutralinos, the inner plot (right panel) shows 
that the parameter /x comes out to be larger than Mz and /1 2 — M 2 2 > M\ hold for most 
of the parameter space. These hierarchies imply that the lightest neutralino Xi is mainly a 
B-'mo and has mass m^o « M\ , while X2 an d the lightest chargino Xi are almost degenerate 
and are mainly PU-inos with mass ~ M2 [30]. Therefore, both these masses do not exhibit a 
significant dependence on Mt, as seen in Fig. [8] (right panel). The heaviest chargino xf an d 
neutralinos X3 4 are mostly higgsinos with mass set by the \x parameter, increasing therefore 
with Mt- 

Notice that, increasing the value of Bt, all the sparticles become linearly heavier since 
they scale as B T - 

Finally, also the Higgs boson masses can be predicted in our scenario. As already men- 
tioned in Section I6~T1 the mass of the lightest CP-even Higgs boson (rrih ) has been computed 
by including the low-energy radiative corrections. In the allowed parameter space of Figl5l 
rrih turns out to be in the range 110 — 120 GeV, thus being testable in the near future at the 
LHC (mainly through the Higgs decay into 2 photons |55j). The heaviest CP-even (H), CP- 
odd (A) and charged Higgs bosons are much heavier. At tree- level, = \iBnj sin 2(3 and 
for m A ^> Mz all these states are almost degenerate, m# ~ m H ± ~ tua- For B T = 20 TeV, 
Fig. [8] shows that tua, mH, m H ± ~ 400 — 470 GeV. For larger values of Bt the masses of such 
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Mt (GeV) M T (GeV) 

Figure 9: Left panel: SUSY muon anomalous magnetic moment contribution Sa^ as a function of 
Mt (taken with its corresponding A along the lines of constant BR(/i — > ey) = 1.2 x 10~ n , 10~ 14 ) for 
S13 = 0,0.2 and Bt = 20 TeV. Right panel: curves of fj, (lower plot) and tan/3 (upper plot), for the 
same inputs considered in the left-panel. 



non-standard Higgs bosons increase, while mu increases by a few GeV due to the logarithmic 
sensitivity to Bt- The Higgs sector is therefore characterized by a decoupling regime with a 
light SM-like Higgs boson (h) and the three heavy states (H, A, H ). 

We conclude this section with a comment on the SUSY contribution to the muon anoma- 
lous magnetic moment 5a M . In Fig. [9] we show the Sa^ behavior as a function of Mt in 
correspondence with different values of BR(/i — > ey) and for Bt = 20 TeV. It is well known 

that 5a ^ is induced by the dipole operator whose dominant contributions go as ~ ^^-r— ' • 

soft 

Since the sparticle spectrum does not change significantly with Mt (Fig. |SJ), da^ directly 
reflects the demeanor of tan/3 (cf. left and right panels). The sign of 5a M is the same as the 
one of n and, in the allowed range Mt >2x 10 s GeV, Sa^ <4x 10~ 9 respects the constraint 
(|36|) . As 5a ^ scales as 1/B T , it decreases for larger Bt- 



7.2 LFV: Model Independent Predictions 

We have already described the structure and parameter dependence of the flavour violat- 
ing SSB mass parameters in Section [5j Here, we intend to further investigate the phe- 
nomenological implications by considering several LFV processes (besides \x —>■ &y) such as, 
/i — > eee, r — > /ry, r — > ey, r — > fififi, r — > eee, r — > /iee r — * e/x/i and // — > e conversion in 
4gTi. The present experimental upper bounds and future sensitivities for the BRs of these 
decays are collected in Table [0 

Let us briefly recall some points related to the computation of such processes. The 
radiative decays £j — > li + y are induced by the effective dipole operator: 

errij {iD$ i o p °j° + iDf^a^ij + h.c.)F pa , (41) 
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where e and F pa are the electric charge and the electromagnetic field strength, respectively. 
The corresponding branching ratios are: 

487r 3 a 

BK(/- 



4 + 7) 



Lr j? 



i 31 



L \2 



+ \D H - 



R\2 



BR(£j ^Pi) 



(42) 



where a = e 2 /(47r), Gp is the Fermi constant, BR(/i 



1, BR(r — > fiVrVfj 



17% 



and BR(r — > ev T v e ) pa 18%. Since LFV resides in the parameters m| and A e , the coefficient 
D L is the dominant one. In the mass-insertion approximation its parametric dependence is: 
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(43) 



where the coupling g can be either $2 o r <?' = y 3/5gi- In the MSSM framework, the dipole 
coefficient D L receives contributions from three different diagram topologies which, under 
certain features of the sparticle spectrum, could cancel against each other. In the latter 
case, these decays would be strongly suppressed even if the sparticle spectrum were not too 
heavy [SB]. This can be realized if, for instance, the A-terms A e are large, or the gaugino 
mass parameters M\ and M2 have opposite sign. In our scenario such peculiar situations 
do not occur, so the suppression of the dipole operators can only arise from large mass 
parameters (i.e. by increasing Bt) or from cancelations inherent to the m| off-diagonal 
entries, as discussed in Section O 

The strict predictions regarding the ratios (133!) and (!34j) can be translated into predictions 
for the ratios of the BRs [9] once Eqs. (142]) and (1431) are taken into account: 
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(44) 



where the results on the r.h.s hold for all the three types of neutrino mass spectrum and the 
numbers in parenthesis apply for IH case, whenever this is different from the others. 

Regarding the 3-body decays tj — > ti&jjtk and ji — > e conversion in Ti, these receive 
contributions from 7- (dipole D L ' R and monopole C L)R operators), Z- (monopole A L,R 
operator) and Higgs- (monopole A L,R operator) exchange diagrams, as well as from box 
(B LL,RR , B LR ' RL ) diagrams. The branching ratios are given by: 
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and the coefficients F NP (N, P = L,R) are 
see e.g. [56J). We have found some discrepancies 



where a k = 1(2), b k = 3(11/4) fork^i (k 
combinations of A N , B NP , C N and A^ 
between our analytical computations concerning the coefficients A L,R with those published 
by some authors [IT] . In Appendix 2 we have reported and discussed our results for A L,R . 
Moreover, in Appendix 3 we present the results for the box coefficients B MN relevant for the 
decays r — ► fiee and r — > e/i/z, calculated to all order in the electroweak-breaking effects. 
For the sake of brevity, we refer to the existing literature (e.g. [UJ and [56] ) for the explicit 
formulas of all the other coefficients, including the one relevant for the CR(/i — > e;Ti). 
Like for the dipole D L , the presence of LFV in the left-handed sector implies that only the 
coefficients A L , B LL , B LR , C L , A L are important. We just recall their parametric dependence: 
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From what has been commented above about suppressing D L in our scenario, and by com- 



paring Eq. (T431) with (f4"6"|) . one can realize that, whenever the coefficients D L get suppressed, 
A L , B LL ' LR and C L undergo the same fate. Regarding A L , one has that A^ — > only if 



ml 



* L )ji — > since these coefficients are insensitive to an overall mass scale increasing. Con- 
sequently, the D L contributions to the 3-body BRs [CR(/i — > e)] are dominant with respect 
to those from the remaining operators, due to the tan (3 and phase-space logarithmic-factor 
[tan/3] enhancement. 

Although the Higgs-exchange diagram contribution also benefits from a tan (3 enhance- 
ment, its numerical relevance requires the Higgs bosons A and H to be significantly lighter 
than the sleptons, charginos and neutralinos EH EU [59]. As already discussed in Sec- 
tion 17.11 in our scenario this does not occur, thus the Higgs-mediated contributions come 
out to be subleading. 

In the aforementioned dipole- dominance situation, the LFV processes under consideration 
can be directly compared with the radiative decays: 
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is the muon capture width 



4.5 x 10 s is the total muon decay width. In addition to these 
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Table 2: Expectations for the various LFV processes from Eqs. (I44|) , f|47l) and (|48l) , assuming BR(/i — » 
e*y) = 1.2 x 10~ n . The results in parenthesis apply to the case of the IH neutrino spectrum, whenever 
these are different from those obtained for NH and QD. 



correlations, we also find that the 3-body decays can be related to each other by using the 
ratios (1511 as follows: 
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(48) 



where the parenthesis enclose the results for the IH spectrum (when different from the 
other cases). In Table [2] we present a synoptic view of the correlation pattern predicted in 
our context, assuming that the present bound on \i — > e r y (Table CD) is saturated, choosing 
s i3 = (0.2) and setting all the remaining neutrino parameters at their best fit points 
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7.3 LFV processes: Numerical analysis 



We are now in position to analyse our numerical results regarding the predictions on the LFV 
decay branching ratios. At this point, we are interested in knowing how and to what extent 
the LFV processes can probe the parameter space of our framework. In other words, will the 
upcoming experimental sensitivities be enough to test the allowed parameter space of Figs. 
and [HP Instead of plotting the contours relative to the various BRs in the allowed (A, Mr) 
space, we have fixed Bt = 20 TeV, selected some values of Mr and performed the analysis 
along one direction, which can be either S13 (and the phase 5) or A. This is representative 
enough and allows us to properly display the main features. 

Consider the effective size of LFV and QFV, which can be parameterized by the following 
dimensionless parameters [HI [8] : 

r |(m~),-,-| , |(mj);il 

* =^#. (i,j = e,ii,T) , 4 = ^fM (i,j = d,s,b), (49) 

L d c 

where m| and m| are the average L and d c masses, respectively. The parameters 8 L,d are 
independent of Bt, while their overall size is determined by the ratio Mr/A. In Fig. [TOl we 
have plotted 8^ and Sfj as a function of s 13 for two selected points of the parameter space 
shown in Fig. EJ (A, Mr) = (4.8 x 10~ 5 , 10 9 GeV) (Fig. M, upper panels) and (A, Mr) = 
(2.4 x 10~ 5 , 10 9 GeV) (Fig. (TUl lower panel). The first point lies in the allowed portion of 
the parameter space studied in Fig. [5] for S13 = (very close to the region delimited by the 
present \i — > e 7 bound), whereas it is excluded for s i3 = 0.2. The second point falls into the 
region excluded by the /i — > e 7 bound, for the NH (and QD) spectrum, irrespective of S13 
(but it is allowed in the IH case for some values of S13). 

The behavior of 8 L,d just reflects that of the relevant FV structure V(m^') 2 V"'" in Eq. (|26|) . 
For example, notice that they scale as 1/A 2 outside the cancelation 'dip'. This figure clearly 
shows that 8^ , 8 d s are insensitive either to si 3 or to the type of neutrino spectrum (cf. the 
upper panels), while 8^ e (8 d d ) gets exchanged with 8^ e (8 d d ) when passing from the NH to the 
IH case. This latter feature is due to the fact that 623 = tt/4 and so the flavours /1 and r 
(or b and s) are indistinguishable. The relative ratios 8^ /8^,, 8 d s /8 d d and 8^ e /8^ e , 8 d d /8 d d 
just reproduce the absolute values of -R23/12 and -R13/12 in Eq. (|34"|) . respectively. As we 
have already deduced from (T26]) . this constant-ratio rule is violated for S13 ~ 0.02, where 
8^ e , 8 d d and 8^ e , 8 d d are strongly suppressed for the NH and IH spectrum, respectively. All 
the above peculiarities are common to the related curves of BR(£j — > ^7), plotted on the 
right of each panel. In particular, the ratios fj4*4l) are remarkably reproduced, except for 
S13 ~ 0.02 where BR(r —> ej) and BR(/i — > ej) undergo a sharp suppression for the NH 
and IH spectrum, respectively. Hence, the point with A = 4.8 x 10~ 5 and S13 ~ 0.02 is 
excluded by the BR(/z — > ej) bound in the NH case, whereas it is allowed in the IH with 
BR(r -»• /17) ~ 3 x 10~ 9 , BR(r ei) ~ O(10^ n ). The point with A = 2.4 x 10~ 5 (lower 
panel) is phenomenologically viable if the neutrinos have IH masses and S13 ~ 0.02. In such 
a case, BR(/x -> ery) < 10 -11 , BR(r -> /x 7 ) ~ 5 x 10~ 8 and BR(r -> ery) ~ O(10- 10 ). 
This is an example where BR(r — > fij) is close to the present bound and /i — > e 7 might be 
unobservable. 

In the examples considered above, the size of the 8 d, s in each family sector is smaller by 
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Figure 10: FV parameters Sfj,(i,j = e,n,r), 5fj,(i,j = d,s,b) and BR(£j — ► ^7) as a function of 
Sx3 for the NH (left upper-panel) and IH (right upper-panel and lower panel). The parameters Mr, 
St and A are fixed as indicated on the top of each panel. For each radiative decay, the present upper 
bound (horizontal solid-lines) and future experimental sensitivity (horizontal dashed-lines) for the BRs 
are shown. 
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Figure 11: Contours of the ratios BR(r — > /x7)/BR(/i — > e7) (dashed lines) and 
BR(r — > e7)/BR(/i — > ej) (solid lines) in the (S13, 5)-plane, for the NH (left panel) and IH 
(right panel) neutrino mass spectra. 



a factor of ~ 3.5 than the one of the corresponding 5 L, s. This results from different com- 
pensating effects. Namely, the squark masses are about three times larger than the slepton 
ones (due to the gluino- induced renormalization effect; see also Fig. |5j), while the (nr-,)^- 
are larger than the (m|)y at the messenger scale Mt (because of the major strong coupling 
contribution). In the first example, 8f s ~ 1.5 x 1CT 3 (left panel). This is the maximal value 
attainable by 5f s in the allowed parameter space of Fig. [51 To perceive the phenomenological 
relevance of such 5$ s , we need to confront it with the gluino and squarks masses. Recall that 
for B T = 20 TeV we got tjIq ~ 900 GeV and M 3 ~ 1.3 TeV (see Fig. EJ. Then, <$£ ~ 10" 3 is 
well below the experimental bound posed by the measured BR(6 — >• 57) [60J. The predicted 
value Sg d ~ 5 x 10~ 4 (obtained for S13 = 0.2) also lies below the bound inferred from K° — K° 
mixing [60] . 

Fig. CD] is aimed to address the dependence of the ratios BR(r — > /i7)/BR(/x — > ej) and 
BR(r — > e7)/BR(/i — > ej) on S13 and the CP- violating phase 5 for a NH (left panel) and 
IH spectrum (right panel). The results do not depend on the ratio M^/A as far as either 
the quadratic or the quartic terms dominate in the LFV entries m|. Regarding the ratio 

Li 

BR(r — > //7)/BR(/x — > ej), it slightly decreases (increases) with increasing S13 for 5 < n/2 
(5 > it/2). However, this ratio blows up at S13 ~ 0.02 and 5 w rr (0) for the NH (IH) 
spectrum. The BR(r — > e7)/BR(/i — > ej) is also sensitive to both s 13 and 5. In fact, for a 
given value of Si 3 , the increase of this ratio with 5 is quite modest in most of the S13 range. 
Still, for S13 ~ 0.02, BR(r — > ej)/BK(fi — > ej) goes to zero (infinity) when 5 = (it) for 
the NH case. Instead, the opposite occurs when the IH pattern is considered since there is 
an interchange of the roles played by BR(/i — > e-y) and BR(r — > ej) (as already discussed). 
In short, barring the range around S13 ~ 0.02, the predictions (jSJ), (H7j) and (|48p are not 
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substantially altered when the effects of 5 7^ are included. 

Finally, we come to a comparative analysis of radiative decays, 3-body decays and fi — > 
e conversion in Ti considering the three types of neutrino spectrum. Fig. [12] shows the 
dependence of the BRs on A for M T = 10 9 GeV and Si 3 = (0.2) in the upper (lower) 
panels, assuming a NH spectrum {mi = eV). 

For Si3 = there is a cancelation dip only in BR(r — > /ry) (left panel), BR(r — > /zee) and 
BR(r — > 3/i) (right panel) at A ~ 3.5 x 1CT 6 , which lies in the (green) region excluded by 
the negative sparticle search, A < 6 x 1CT 6 (cf. also Fig. [5j upper panel). This dip originates 
from the cancelations of the quadratic and quartic terms in m| discussed in Section |5j In 
BR(/t — > cy) this would take place at a value of A smaller by a factor of ~ (6/78) l l 2 ~ 0.3 
[see Eq. ( 12"S|) ] and so, would fall into the (grey) non-perturbative range. For s 13 = 0.2 the 
dips are instead present for all the BRs, as we expect on the basis of Eq. ( |2"9|) . Outside the 
cancelation regions, the relative ratios of BR(£j — > ^7) are those announced in Eq. (jUJ) and 
the BR(£j — > iiikik) are correlated with BR(£j — > ^7) according to Eq. fj4"7l) . For comparison, 
we have also plotted all these BRs for the case of the QD spectrum with mi = 0.3 eV which 
can be obtained by 'continously' rising the mass mi in the NH case beyond (Am 2 ^) 1 / 2 [cf. 
Eq. (j!9p ] . In the QD case the non-perturbative range extends much above the one relative 
to the NH (grey) so that the perturbativity lower bound on A (indicated by the vertical solid 
line) is larger by about one order of mag nitude, A > 3 x 10~ 5 . Notice that a more restrictive 
lower bound is imposed by sparticle searches (green), A > 5 x 10~ 5 . The cancelations occur 
at approximately the same A (inside the excluded regions) for all the BRs. For A > 5 x 10~ 5 , 
the curves corresponding to the NH and QD are superimposed and so the two scenarios are 
not distinguishable [cf. Eqs. (|25|) and (|26|) ]. 

Finally, Fig. [U reveals that, for M T = 10 9 GeV and assuming a NH spectrum with 
S13 ~ 0, only BR(/i — > ej), BR(r — > /17), BR(/i — > 3e) and CR(/i — > e; Ti) are within the 
future experimental sensitivities for 5 x 10~ 5 <A<7x 10 -5 , while if A > 7 x 10~ 5 only 
ji — > e7 and CR(/i — > e; Ti) could be accessible. All the other LFV processes would be 
undetectable in the allowed A range. Then, if for example the MEG experiment detects 
BR(/i -> ery) at the level of 8 x 10" 12 then BR(r ^7), BR(/i 3e) and CR(/i -> e; Ti) are 
expected to be ~ 2 x 10 -9 , ~ 5 x 10 -14 and ~ 4 x 10 -14 , respectively. The case of the QD 
spectrum is similar: in the range A ~ (6 — 7) x 10 -5 the /ie LFV decays are observable but 
BR(r — > fij) is predicted to be < 10 -9 , and for larger A only fi — > e7 and fi — > e conversion 
would be visible. 

By switching on S13, LFV is enhanced in the /ie and re sectors, but it is not essentially 
altered in the r/i sector. As a consequence, the lower bound on A imposed by the present 
limit on BR(/a — > e^y) is more restrictive (cf. Fig. [5]) and consequently, the BRs of the r/i 
sector are penalised. For instance, detecting /1 — > ej with a BR around 8 x 10~ 12 would imply 
the possibility to measure also BR(/t — > 3e) and CR(/i — > e; Ti) at the level of 5 x 10~ 14 and 
4 x 10~ 14 , but BR(r — > /ry) would be ~ 10 -11 (well below the planned future capability). 
This conclusion holds for both the NH and QD spectrum. 

Fig. [13] presents a similar analysis for the IH spectrum with m^ = eV. For comparison, 
we also show the QD case which can be obtained by pushing the mass m% to values larger 
than (Am^) 1 ' 2 thus recovering a case very similar to the QD one reported in Fig. [T2j When 
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Figure 12: Left panels: branching ratios of the radiative decays £j — ► ^7 and [i — ► e conversion in Ti, 
as a function of A. Right panels: branching ratios of the three-body decays £j — > 3£j and r — > //ee , e/x/i. 
All the results are shown as a function of A for Mt = 10 9 GeV and Bt = 20 TeV taking S13 = 
(upper plots) and S13 = 0.2 (lower plots). The left (right) vertical line indicates the lower bound on A 
imposed by requiring perturbativity of the Yukawa couplings Yt,s,z when mi = 0(0.3) eV [NH (QD) 
neutrino mass spectrum]. The regions in green (grey) are excluded by the rug > 100 GeV constraint 
(perturbativity requirement when mi = OeV). 
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Figure 14: As in Fig. Q2 but for M T = 10 12 GeV. 
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Figure 15: As in Fig. O but for M T = 10 12 GeV. 



•5i3 = (upper panels) the main difference with the NH case is the fact that all BRs exhibit 
the suppression dip in the perturbative range (in the r/x sector this takes place for a smaller 
A with respect to the /xe one). However, due to the choice of Bt these dips fall into the 
(green) range excluded by sparticle searches. Outside these cancelation regions the values 
of all these BRs (and the related correlations) are essentially the same as those obtained for 
the NH case leading to the same phenomenological implications. 

The IH scenario offers the possibility to reverse the dominance of the /xe LFV over the r/x 
LFV. Suppose that the effective SUSY breaking is increased up to B T > 70 TeV so that all the 
sparticle masses are enhanced by a factor of 3.5. Then the BR(/x — > ej) dip would lie in the 
allowed A rangd^and BR(r — > fij) would be brought close to the present bound. As a result, 

16 The effect of increasing Bt does not modify the perturbativity bounds, but for Bt <D 50 TeV the 
constraint from the negative sparticle searches disappears, as the comparison between Fig. [5] and [6] has 
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fi — > ej, fi — > 3e and fi — > e conversion might be invisible, whereas BR(r — > /ry) ~ 5 x 10~ 8 
could be measured [perhaps together with BR(r — > fiee) ~ 5 x 1CT 10 ]. 

By switching on S13 = 0.2 (lower panels) we find, as for the NH case, that all the 
LFV processes undergo a dramatic suppression at the same A which, however, occurs in 
the excluded range. Outside this interval, the relative ratios of the BRs closely follow the 
predictions^ OS]), (jSJ) and (T4"8"j) . In such a case, the planned experimental sensitivities would 
allow to test only fi — » ej, fi — ► 3e and — ► e conversion. 

Figs. [14] and [15] contain the results of a similar analysis performed by taking My = 
10 12 GeV for the NH and IH spectrum, respectively. The increasing of Mt implies an overall 
decreasing of all the BRs due to a heavier sparticle spectrum (cf. Fig. [8]) with respect to 
the case with Mt = 10 9 GeV. We again find that for any A (barring the cancelation ranges) 
the several BRs are correlated according to the model-independent predictions contained in 
Eqs. OS]), (j47J) and (|48j) . The features of the cancelation dips are those discussed in Section [5] 
and already encountered for smaller Mt- Consider the NH case (Fig. [T4"j) . At S13 = and A 
within (3 — 4) x 10 -2 , fi — > e^y, r — > fij, fi — > 3e and fi — > e conversion are in the reach of the 
future experiments in both the NH and QD cases. From A ~ 4 x 1CT 2 up to 1.5 x 1CT 1 only 
ji — > e'j and fi — > e conversion could be detected. For s 13 = 0.2 and A > 0.1 only fi — > ej, 
— >• 3e and fi — > e conversion are accessible. We arrive at similar conclusions observing 
the case of IH spectrum displayed in Fig. [15) the \xe sector is more favoured with respect to 
the T/i sector, where only the r — > /ry decay could be detected as far as S13 ~ 0. However, 
this conclusion could be contradicted for larger Bt and S13 = 0. Consider B T ^ 50 TeV 
so that the sparticle spectrum increases by a factor of 2.5. Then, in the suppression dips 
of the fie LVF processes (see upper panels in Fig. ITS]) . BR(r — > ftj) would be close to the 
present bound and BR(r — > fiee) ~ 10 -9 . Therefore, the future detection of both r — > fi'j 
and r — > fiee, together with a slepton and squark spectrum in the range 400 — 650 GeV and 
1.8 — 2.5 TeV, respectively, would point to a IH neutrino masses and si 3 ps 0. 

8 Summary and Conclusions 

The future perspectives to detect signals of new physics mostly rely on the observation of 
sparticles at the LHC or LFV decays at e g. the B-factories [38] , the incoming MEG experi- 
ment [36], the Super Flavour factory [61] or the PRISM/PRIME experiment at J-PARC [46]. 
It is, therefore, extremely important to motivate and suggest theoretical scenarios which can 
be tested in more than one direction. In this paper we have presented and discussed in detail 
a supersymmetric 577(5) version of the triplet seesaw mechanism in which the triplet are 
messengers of both L and SUSY breaking. The key-points of our model can be outlined as 
follows: 

• The tree-level exchange of the triplets T generates neutrino masses, so flavour violation 
is induced also in the lepton sector of the SM, as required by the observation of neutrino 
oscillations. All the LFV effects are parameterized by a single flavour structure Yy; 

already demonstrated. 

17 The fact that in Fig. [13] BR(r — * 3/i) is slightly larger than BR(/i — > 3e), in contradiction with the 
predicitons reported for example in Table [21 is due the Y" T -driven RG running of the neutrino mass m„ which 
concerns mainly the entries ri [i = e, /1, r), and is more sizeable for S13 7^ with IH spectrum. 
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• The quantum-level exchange of the 15-states T, S and Z generates all the SSB mass 
parameters of the MSSM via gauge and Yukawa interactions. Their mass scale is 
determined by the effective SUSY breaking scale Bt- Flavour violation is induced 
in the mass matrices m| and m|, and in the trilinear terms A e and A^, by the 
Yukawa couplings Y T (and, according to the 577(5) relation f[T2"j) . also by Y s and Y z )- 
Therefore, the flavour structure of the SSB parameters is provided by Yy or, according 
to Eq. ffTSl) , by the low-energy neutrino parameters. 

• The number of free parameters is three: the triplet mass Mr, the effective SUSY 
breaking scale Bt and the unflavoured coupling constant A. 

These aspects make the present scenario highly predictive since it relates neutrino masses 
and mixing, sparticle and Higgs spectra, lepton and quark flavour violation in the sfermion 
masses and electroweak symmetry breaking. We have performed a complete analysis of the 
parameter space spanned by Mt, Bt and A taking into account the present experimental 
constraints on the above physical observables (see Figs. [5] and [6]). This has demonstrated that 
there is a region in the parameter space where our framework is compatible with experiment. 
In particular, despite the very constrained structure of our SSB mass parameters, EWSB can 
be radiatively realized analogously to the conventional MSSM case. Our predictions allow 
us to further test the allowed parameter space as follows: 

• Regarding the MSSM sparticle spectrum, we predict that the gluino is the heaviest 
sparticle while, in most of the parameter space, i\ is the lightest. However, the gravitino 
is lighter than the MSSM sparticles. For instance, if Bt = 20 TeV, the squark and 
slepton masses lie in the ranges 800 — 950 GeV and 100 — 300 GeV, respectively. 
The chargino masses are m^± ~ 320 GeV and m-± ~ 450 — 550 GeV. Moreover, 
m^o ~ 190 GeV, m^o « m~± and m^o »s m~±. These mass ranges are within the 
discovery reach of the LHC. Increasing the parameter Bt implies a linearly heavier 
spectrum. The measurement of only a few sparticle masses will provide a hint on the 
value of the effective SUSY breaking scale Bt and a test of the correlation pattern 
shown in Fig. [HJ 

The mass range of the electroweak sparticles implies that the SUSY contribution to 
the muon anomalous magnetic moment Sa^ never exceeds the maximum value of the 
interval fl36|) . Moreover, since Sa^ > 0, the discrepancy between the experimental 
determination and SM prediction is alleviated in our model. 

• The Higgs sector is characterized by a decoupling regime with a light SM-like Higgs 
boson (h) and the three heavy states (H, A and H ± ) with mass rriH,A,H± ~ 450 — 
550 GeV (again, for Bt = 20 TeV). These masses increase almost linearly with Bt- 

• We have considered several LFV processes: fi — > eX, // — > e conversion in nuclei, 
t — > eY and r — > jiY (X = 7, ee, Y = j,ee,nn). Our framework is characterized 
by peculiar LFV and QFV patterns [intimately related to each other; see Eq. fl33l ]. 
which are mostly determined by low-energy neutrino masses and mixing. The size of 
QFV, when confronted with the coloured sparticle spectrum, is well below the present 
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phenomenological bounds extracted from b — > 37 and X° — K° mixing, etc. Therefore, 
QFV processes do not constrain our parameter space. 

Concerning LFV, we stress that strict predictions have been obtained for the relative 
branching ratios of the radiative £j — > la [see Eq. ( 1441 ] and 3-body £j — > Idk^k decays [see 
Eq. (l48l) ]. The latter, as well as \i — > e conversion in Ti, are also correlated with the radiative 
decays as shown in Eq. (|47[) . All these results are model-independent in the sense that they 
do not depend on My, Bt or A (they are given only in terms of the low-energy neutrino 
observables) . If the present bound on BR(/i — > cy) is saturated, the branching ratios of the 
remaining LFV processes are predicted as shown in Table El where the three types of neu- 
trino spectrum have been considered. The experimental signatures of our scenario crucially 
depend on the value of the lepton mixing angle #13: 

Tiny 6 13 : The analysis has shown that, in the allowed parameter space, the future exper- 
imental sensitivity will allow to measure at most BR(ii — > ej), BR(/i — > 3e), BR(r — > fj,j) 
and CR(/i — > e Ti) according to the relations (Pill) and (PJ7]). In particular, being BR(r — > 
/i7)/BR(/x — > ej) ~ 300, BR(r — > 7x7) is expected not to exceed 3 x 10~ 9 , irrespective of the 
type of neutrino spectrum. All the decays r — > iii^k would have BR < O(10~ u ). 

Sizeable 1S : If S13 is close to the upper bound fl23|) . the r/i sector is hardly accessible and 
only the decays fi —>■ &y, /J — > 3e and fi — > e conversion in Ti can be observed in the future. 
This conclusion holds for the NH, IH and QD neutrino spectra. For instance, if the MEG 
experiment measures BR(/i — > cy) ~ 8 x 10 -12 , then BR(/i — > 3e) and CR(/i — > e Ti) are 
expected to be ~ 5 x 10" 14 and ~ 4 x 10" 14 . The latter is in the reach of the PRISM/PRIME 
sensitivity. Values of S13 ~ O(10 _1 ) will be explored soon by several neutrino experiments 
like MINOS, OPERA and Double Chooz g2]. 

Deviations to this very specific model-independent pattern occur when (nx|)y cancel (see 
discussion in Section EJ). This can be the case if S13 ~ 0.02(— 0.02) and the neutrino spectrum 
is IH (NH). Then, all the /ie LFV processes might be invisible and r — > /J7 detected with 
a BR close to the present bound (taking Bt = 20 TeV). Moreover, r — > /xee would be 
O(10~ 9 — 10~ 10 ) and the remaining r decays below 10~ 10 . Notice that, increasing Bt the 
BRs are suppressed since they scale as 1/B T . This specific value of S13 is in the sensitivity 
range of Neutrino Factories [53] . 

Alternatively, (m|) Me could be vanishing because of a cancelation between the quadratic 
and quartic Yukawa contributions [see Eq. (|30|) ]. For the IH case, if Bt ^ 50 TeV and S13 
is very small, all the /ie LFV processes would be strongly suppressed whereas the t\l sector 
would be favoured with BR(r — ■> fij) ~ 5 x 10 -8 and BR(r — > fiee) ~ xlO -9 . This scenario 
is correlated with a slepton and squark spectrum in the ranges 400-650 GeV (in the limit of 
the LHC detection capabilities) and above ~ 2 TeV (within the LHC reach), respectively. 

Given the increasing interest on the problem of finding a possible relation between lep- 
togenesis [61] and low-energy neutrino physics [65], we would like to comment on this issue 
in the framework of our work. Within the present version of the supersymmetric triplet 
seesaw mechanism, leptogenesis can be realized by considering that the soft bilinear term 
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Bt produces a mass splitting between T and T, leading to resonant leptogenesis [66]. For 
this to work, Bt must be around the electroweak scale. In our case, since Bt > 20 TeV, the 
BAU turns out to be too small. Still, leptogenesis could be made effective either by adding 
an additional pair of triplets [H] or by including heavy singlet neutrinos [67]. The former 
case would imply the appearance of one more flavour source Y^ which, depending on its 
size, could have some impact in the predictions of our framework. Instead, we would like to 
comment on the second possibility which requires heavy neutrino singlets (with mass M^) 
coupled to the lepton doublets through the Yukawa couplings Y^. To maintain the predic- 
tions made along this work we must require that the singlet contribution to neutrino masses 
is much smaller then the one generated by the triplet i.e., Y 2 N <C XYtMn /Mt- Moreover, 
Yat Yt is also required to suppress LFV arising in the SSB parameters from the sin- 
glet exchange at the quantum level 18 ! (notice that the neutrino singlets could couple to the 
spurion field X) . For the purpose of leptogenesis, two different situations can be envisaged: 
M/v > Mt or < Mt- In the former case, the CP-asymmetry generated through the 
decay of the triplets into two leptons is directly proportional to Y 2 N /M N [67], so it turns out 
to be very tiny (Y 2 N <C XY T M N / M T )- On the contrary, if M N < M T the CP- asymmetry is 
weakly sensitive to Y N [67], therefore a viable value for the BAU can be achieved. 

We conclude our discussion by remarking that our scenario is not only extremely pre- 
dictive but it can also be tested in view of the present and future programmes of LFV and 
neutrino oscillation experiments. 
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Appendix A-Extracting the SSB terms from wave func- 
tion renormalization 

In this Appendix we derive the general expressions for the soft supersymmetry breaking 
scalar masses, the bilinear and trilinear couplings at the messenger scale M. We employ a 
generalization of the method suggested in Ref. [30] and subsequently presented in Ref. [31 J . 
Consider the case in which the scales of SUSY breaking and its mediation to the observable 
sector are determined by the vev of the auxiliary and scalar components of a chiral singlet 
superfield X, (X) = M + 6 2 F (in the following M is taken to be real and |P| < M 2 , for 
consistency of the method). The leading contributions (at lowest order in P) to the SSB 
terms arise from X-dependent wave-function renormalizations of the chiral superfields Zq. 

18 For a discussion related to this in the context of EWSB see [55] . 
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The effective Lagrangian readj^l 

£ = J d*dQtZ Q (X,X^) a b Q b +\J d 2 9(/i ab Q a Q b + f abc Q a Q b Q c ) +h.c. 



(A-l) 



where fj, a b and f a t, c are superpotential mass parameters and dimensionless coupling constants, 
respectively. The wave-function renormalization Zq is a hermitian matrix which depends on 
\X\ = \fXJO. Its ^-expansion at X = (X) is given by 
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After expressing Eq. (lA-lj) in terms of the canonically normalized superfields Q' as 

e =( 1+ ¥if 92F ) z « /2s '' (A - 3) 

we can extract the SSB masses for the scalar component of Q from the quartic terms 9 2 9 2 
in the first integral and the SSB bilinear and trilinear couplings from the quadratic terms 6 2 
in the second and third integral, respectively 
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In order to find the explicit expressions for the SSB parameters at an energy scale /i we recall 
that the /x-dependence of Zq is expressed by the RG equation: 
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where t = ln/i and 7q is the matrix of anomalous dimension. By defining Zq 2 = 1 + SZq 2 , 
where SZq 2 encodes the quantum corrections, Eq. (1A-5|) reads 
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(A-6) 



19 Hcre the index a labels either the superfield and its associated 'charges' or only the 'charges'; the context 
should make clear the case. 
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at the lowest order. By following the procedure outlined in Ref. [31] , we have to formally 
integrate the above equation to obtain SZq in terms of the anomalous dimension. After- 
wards, the solution can be plugged into the expressions (1A-4|) to extract the SSB terms. 
For the sake of brevity, we do not report all the intermediate steps (which can be easily 
performed) and, instead, give the final expressions at n = M: 
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where 7g(M) [7g(M)] is the anomalous dimension above (below) the mass scale M, A 7 g(M) 
7g(M) — 7g(M) and (A^g)7g(M) means considering the difference of the beta-functions 
of the couplings contained in 7g(M) above and below M. Notice that our result f]A-7|) for 
differs from the one obtained in Ref. [3T] because an extra term proportional to the 
commutator [7q,7q] appears in that work, which is manifestly inconsistent for a hermitian 
quantity such as m|>. We believe that the appearance of such a commutator is due to the 
improper definition of the RG equation for Zq given in Eq. (3.5) of Ref. [3T] . 

In the following, we provide the explicit expressions for the relevant anomalous dimen- 
sions needed to extract the SSB parameters from Eq. (I A- 71) in our specific framework. The 
anomalous dimensions below the scale M are: 
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Regarding the expressions for the RG equations at one loop in the MSSM framework with the 
(15 + 15) SU(5) representation we refer to Ref. [9]. Finally, we obtain the explicit formulas 
given in Eqs. ffTBl) and ffTTl) . 



Appendix B-Coefncients of the ijiiZ operators 

In this Appendix we compute the coefficients A L)R of the monopole operators 

gz(AfAa% + A%i\an) + h.c.)^ , A$*> = A^ R) ' C + A$ R) > n (B-l) 

where gz = g2/c\v ( c w — cos9w, Q\v is the weak mixing angle) and A L ^ R ' ,C [A L ( RS) ' n ] stand for 
the contributions from the chargino/sneutrino [neutralino/charged-slepton] loop diagrams. 
The two-component spinor notation is used such that, for example, &i(£f) is the left-handed 
(right-handed) component of the lepton i field (i = e,{i,r). Different one- loop results for 
such coefficients have been presented in the literature. For instance, the authors of Ref. [TT~] 
provided an all-order calculation in the electroweak breaking effects, while the authors of 
Ref. [56J performed a lowest-order calculation. We found dramatic numerical discrepancies 
between the two aforementioned results, which cannot be ascribed to the approximation used 
in Ref. [56] . Moreover, the authors of Ref. [59] have recently re-evaluated the contributions 
to A L ' R and claimed to have found additional contributions disregarded in Ref. [TT] . We have 
independently performed the all-order computation to compare with the previous results and 
to clarify this issue. The notation of Ref. [TT] has been adopted to define the mass eigenstates 
of the charged sleptons £ x , (X = 1, . . . 6) , sneutrinos u x , [X = 1, ... 3), charginos x A , (A = 
1,2) and neutralinos x° A , (A = 1,...4). The corresponding (unitary) mixing matrices are 
denoted by U Xi , U Xi , (Ol)a«, (Or)a« and (Otv)a/3 (where a{(5) is the current-basis index for 
charged or neutral gauginos/higgsinos). The relevant interactions of sleptons/leptons with 
charginos and neutralinos are: 

C = liCf Ax v x l\ + e\Cl Ax u x x A + hN^JxXA + <N t L Ax i x XA + h.c, (B-2) 



where 



c ux - -92{0 R ) A1 U Xi , C^ AX - g 2 —7= -{0 L ) A2 U Xi , 

V 2,m w cos p 



N R AX = -»L \ -[{0 N ) A2 + {0 N ) AL tmB w ]U e Xi + - ^ (0 N ) A3 U x(i+3) , 




N iAX = Jk 2 ( N)Aitan6 w U x(i+3) + — — — — ~(0 N ) A3 U x{i+3) 



(B-3) 



m w cos (3 ' 
mw cos (3 ' 

The interactions of charginos and neutralinos with the Z boson are the following: 

£ = ~gz (WA^xtRAB - Ta^XbLab + Fa^X° B N ab ) Z„, (B-4) 

where 



Ral 



c 2 e w (0* R ) A1 (0 R ) m + (i - sj w ) (0* R ) A2 (0 R ) B2 



Lab — R*ab\r-^l, 



Nab = \ l(0* N ) A3 (0 N ) B3 - (0* N ) M (0 N ) m \ (B-5) 
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The chargino contributions are: 



A L,c 



^iBX^jAX 
16^2 



F{%X A i> x , x xbi>x) + m MT 

E>* A 
Ka 



AB ' 



in 



SabI^- Z\ 



H(x* Ai>Y ) + \n-Vk 



aR,c 



A L ' C \ , 
A ji \l* 



.R- 



;b-6) 



;b-7) 



here A,B = 1, 2, X 



1, 2, 3 and Z[ 



Tf 



2 , (Tl 



1 T^3 

~2> J ^fl 



(A summation over repeated indices is understood.) In the above equations the first and 
second terms come from the diagrams in which the Z boson line is attached to the chargino 
line, the third one where it is attached to the sneutrino and the fourth term comes from 
the wave function renormalization. For completeness, we have also displayed the terms 
proportional to ln(// 2 /M|) in Eq. (1B-6I) . coming from the divergent diagrams, where /1 is the 
renormalization scale. Obviously, such terms cancel out. As for the argument of the loop 
functions we have adopted the convention x afe = m^/m 2 , then e.g. x XA u x = M? A /ml x . The 
loop functions are defined as follows: 



F(x,y) 
G(x,y) 
H(x) 




Using the relation H(x) = F(x,x) — 2xG(x,x), one can easily verify the validity of the 
Ward-Takahashi (WT) identity in the SU(2)w x U(l)y unbroken phase, which entails the 
vanishing of the coefficient, A L( - R ^ c \ Vl=V2=0 = 0. By exploiting this identity in Eq. (1B-6I) . the 
above expressions (1B-6|) simplify to: 



.4 



L,c 
ji 



fiR riR* 
^iBX^jAX 

167T 2 



m \ (n*\ ^{ x xavxi x xbvx) 
[y r) A2\y r) B2- 



-(0* l ) A2 (O l )b2 — * A 2 —G(x% A v 



.1 \ 



in 



XA^xi ^Xb^x 



VX 



A 



R,c 



-A z 



l\l^r- 



(B-9) 
(B-10) 



The result obtained for A L in (1B-9[) coincides with that of Ref. [Hj and is consistent with 
the one in Ref. [56] • The formulas (1B-6| IB- 71) are also in agreement with those reported 1 ' 



20 In fact, the agreement between Eqs. (|B-61 IB-7|) and the corresponding formulas in Ref. [59] does not 
regard the constant terms in the loop-functions F(x, y) and H(x). Nevertheless, such terms do not contribute 
because of the unitarity relations. 
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in Ref. (59], which, however, have not been reduced to the form flB-9t IB-lOj) . We observe 
that A% c < Al ,c because of the Yukawa coupling suppression (this coefficient has been set 
directly to zero in [TT]). 

The coefficients A™ R ' ,n from the neutralino-exchange contributions are given by: 



A L,n ^iBX^jAY i , 
A J = — ~ir < Oxy 



16vr 2 
Sab 



n* ab - 



F{x 



,x- 



+ ln- 



\ \t x • " xd'x ' "' \ri » r M^ A My, B 



2 r 
S\v°XY 



H(x 



XA"X . 



+ In 



Ml 



— n e * n e - 

u Xk u Yk 



+ N 



AB 2 

m 



G{x 



XA?X ' X XB~tx 



+ In 



H(x- ) + ln-&, 
-tOabOxy^l o 



A R, n 
A 0i 



Nt BX N AX j , 

J Oxy 



16tt 2 



F(x 



N 



+ ln- 



AB ' 



+ N 



AB 2 

mi x 



G ( x XAe x > x XBlx) 



'AB 



s 2 w 5 X Y (h(x^ x )+\jij^ 



H(x- ) + ln^ 

e V XAlx> ]_ 



TTl* TT 

U yuU- 



ln^ 



where A, B — 1, . . . 4, X = 1, . . . 6, k — e, r and the loop function /(x, ?/) = 1+F(a;, y). The 
first and second terms derive from the contributions with the Z attached to the neutralino 
line, the third and fourth from those with the Z attached to the slepton line, and the fifth 
one from the wave-function renormalization diagram. By using again the WT the above 
expressions simplify aJ^l: 



.4 



L,n 



N&xNfo 



16tt 2 



>XY 



N * F ( X XAix> X XBlx) i N ^XA^XB 

N AB ~ V ISj- 



<AB~ 



m 



G{x 



XAl-x ' X XbIx 



+$AB 



ji* Tjl 

>Xk u Yk 



I{x 



x ixi v ) 



xaIy ' ^i x e Y 



H(x 



A 



R,n 
ji 



NlhxN^Y j s 



16n 2 



XY 



4 

Nab 



>XY- 



X.aIx' 



4 



(B-12) 



F ( x xaIx> x ZbIx) , at* M xaM Xb 



m 



G(x 



XAix ' X XB~tx 



21 Using the simplified formulas (|B-91 IB-10)) and (|B-121 IB-13P is more convenient also because cancelations 
are already accounted for. Needless to say that the constant numerical addenda appearing in the loop 
functions (IB-8j) do not contribute to the final amplitudes because of unitarity of the mixing matrices U v and 
U e . Still, they are essential to prove the WT identities and then to yield the simplified formulas (|B-9l IB-10|) 
and (|B^T2l |B^13]) . 
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-o AB u Xk u Yk > . 

(B-13) 

Our results flB-111) are compatible with those of Ref. [59] (see, however, the comment in 
Footnote [20]), and are also compatible with Ref. [56J. Instead, our final expressions flB-121) 
differ from those of Ref. [llj . because of the third and fourth (third) terms in A L,n (A R ' n ). 



Appendix C-Box coefficients for the lj — > lilk^k (i 7^ k) 
amplitudes 

We consider the four-fermion operators which are relevant for the amplitude of the LFV 
decays £j — > Idk^k'- 



,) (B% k £ k a,i k + Bf* £ c k aJt) + h.c. (C-l) 



Each coefficient B MN (M, N = L,R) receive contributions from box diagrams with either 
charginos/sneutrinos or neutralinos/charged-sleptons exchange (see Fig. [T6i) : 



r>MN 



R MN(c) j~)MN(n) 
n ji;k ~r n ji;k 



(C-2) 



We have obtained the following results valid for % ^ k: 



R LL(c) 
n LR{c) 



-J^(M? A ,M^ B ,ml x ,ml Y ) Cf Ay CfXx^kBX^kBY + ^flx^BX^kAY^kBY 

-J<i(M~, A , M^ B ,m 2 - x ,m 2 Y )C R AY CfX x C kBX C kBY 



B 



LL(n) 
ji;k 



in ~ 

XBi v x i Vy 



'2 \(~iK /^iR* (~<L r~iL* 
m v Y '^iBX^'jAX Ly kAY Ly kBY 



B 



LR(n) 
ji;k 



+ 2 M x^xb^(M| a , M 2 B ,m} x ,ml) 



t\tR tvtR,* t\tR t\tR* 
iV iBX 1 ^ jAX ly kBY ly kAY 



+N t R BY N^ x N k H BX N^ Y 
~J 4 (Mf A ,M| B ,mf x ,mf y ) [N&yNfccN&yN&x - N^N^N^yN^y 
+N t R B Y NfX x N k L X Y N k L BX 



T R* 



R 



T R* 



2 m? 



\ P ' ' ' "t-x ' m i Y ^iBxNjAxNkByN kA Y 



n ji;k 
n RR(n) 



r>LL{c) I 
D ji\k \L^R, 



ji;k \L<-tRi 



r)RL(c) 

n ji;k 
j-)RL{n) 
&ji;k 



R LR{c) I 
D jv,k \L^Ri 

n LR{n) 1 
" n ji;k \L^R, 



(C-3) 
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where the coefficients N*^x , Cfj^ have been defined in Eq. (IB-2[) and the loop integrals are 
given as: 



h(m l} . . . ,m 4 ) 



d 4 k 



16tt 4 7 (k 2 - mf) . . . {k 2 - ml) 
Ji(ml,...,ml) = -^7/713 — d \ro — 3. (C-4) 



16vr 4 y (A; 2 - m?) . . . (A; 2 - ml) 

The box coefficients for the decays r — > fiee and r — > e/i/z correspond to the replacements 
(j, i, k) — > (r, /i, e) and (j, i, k) — >■ (r, e, //), respectively. Our results for these coefficients are 
in numerical agreement with those of Ref . [56] . 
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Figure 16: Box diagrams relevant for the computation of the coefficients B^ k (first and second rows), 
Bj^. (third and fourth rows), B^j, (fifth and sixth rows) and B^ L k (seventh and eighth rows) [see 
Eqs. CU) GDI 
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